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Summary 


A survey of studies conducted since 1914 on the use of equivalent-plate 
stiffnesses in modeling the overall, stiffness-critical response of stiffen ed 
plates and shells is presented. Two detailed, comprehensive derivations of 
first-approximation equivalent-plate stiffnesses are also presented that 
are based on the Reissner-Mincllin-type, first-order transverse-shear 
deformation theory for anisotropic plates. First, a derivation based purely 
on static and kinematic equivalence between a stiffened plate and its 
homogenized equivalent is presented, followed by a derivation based on 
equivalence of the strain-energy density. In both derivations, the stiffener 
members are modelled as beams that are shear deformable within and 
transverse to the plane of the plate, consistent with the classical 
continuum mechanics representation of solids. Additionally, each stiffener 
is presumed to be constructed, at most, in a nonhomo geneous manner 
from orthotropic materials with one axis aligned with the stiffener axis 
and the other two axes aligned with the cross-sectional axes. This 
presumption allows the computation of equivalent-plate stiffnesses for 
stiffened panels such as those in which the stiffener caps are reinforced 
with high-strength, pultruded rods. Consistent with a first-approximation 
analysis, inplane bending of the stiffeners and total compatibility between 
the plate skin and the stiffeners are neglected. 

Equivalent-plate stiffness expressions, and a corresponding symbolic 
manipulation computer program, are also presented for several different 
stiffener configurations. These expressions are very general and exhibit 
the full range of anisotropies permitted by the Reissner-Mincllin-type, 
first-order transverse-shear deformation theory for anisotropic plates. 
The expressions presented in the present study were also compared with 
available, previously published results. For the most part, the previously 
published results are for special cases of the general expressions 
presented herein and are almost in complete agreement. Analysis is also 
presented that extends the use of the equivalent-plate stiffness expressions 
to sandwich plates with nonidentical, anisotropic face plates, and 
expressions for equivalent-plate thicknesses are presented. 
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Symbols 


The primary symbols used in the present study are given as follows. Additional symbols are 
defined in Tables 1-9. 


a 


dimension of hexagonal beam grid (see figure 21), in. 


A s area of prismatic stiffener cross-section (see figure 10), in 2 

A cell area of basic cell forming beam grid (see equation (39a) and 

figures 14, 16-18, 21, and 23), in 2 


A n , A 12 , A 16 , A 22 , A 26 , A 66 equivalent-plate membrane stiffnesses (see equation (22a)), 

lb/in. 

A 44 , A 45 , A 55 equivalent-plate transverse-shearing stiffnesses (see equation 

(22b)), lb/in. 


core core 

A n , A,, , A 


core 

16 9 


A 22 , A) 


contribution of sandwich-plate core to equivalent-plate 
membrane stiffnesses (see equation (47a)), lb/in. 


A 


plate 
11 9 


A 


plate 
12 9 


plate 

A , 6 , 


A 


plate 

22 9 


A 


plate 

26 9 


plate 

A 66 


plate contribution to equivalent-plate membrane stiffnesses 
(see equations (23a) and (44a)), lb/in. 


A' u , A ' 12 , A ' 16 , contribution of sandwich face plate to equivalent-plate 

A'T, A'T- aC membrane stiffnesses (see equation (46a)), lb/in. 


A A 

■^44 5 -^45 


contribution of sandwich-plate core to equivalent-plate 
transverse- shearing stiffnesses (see equation (47d)), lb/in. 


A 


plate 
44 ’ 


plate plate 
A 45 , A 55 


plate contribution to equivalent-plate transverse-shearing 
stiffnesses (see equation (23d)), lb/in. 


A' 


plate 
44 9 


A' 


plate 
45 9 


A' 


plate 

55 


contribution of sandwich face plate to equivalent-plate 
transverse-shearing stiffnesses (see equation (46d)), lb/in. 


stiffener 

A 1 1 9 

stiffener 

A 12 ? 

stiffener 

A 16 , 

stiffener 

stiffener 

. stiffener 

A 22 9 

A 26 9 

a 66 


stiffener contribution to equivalent-plate membrane stiffnesses 
(see equation (23a)), lb/in. 


stiffener 
^ 44 9 


stiffener 



A 55 


stiffener contribution to equivalent-plate transverse-shearing 
stiffnesses (see equation (23d)), lb/in. 
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B 


B 1 1, B 12 , B 16 , B 2 2, B 2 6, B 66 


plate plate plate 

^11 9 12 9^16 9 

plate plate plate 

^22 9 D 26 9 D 66 


,p!»e , Pla» 

^11 ’ ^ 12 9 ^ 16 9 

P late plate plate 

^22 t D 26 9 *-> 66 


stiffener stiffener stiffener 

^11 9 AD 12 9^16 9 

stiffener stiffener stiffener 

^22 9 D 26 9 ^66 


[0 


[4] 


[&] 

d s 

Du, D 12 , Di 6 , D 22 , D 26 , D 66 

core core core 

D„ , D 12 , D I6 , 


< ^ core core < ^ core 
U 2 2 9 ^26 9 -L ^66 


plate plate plate 

D D D 

1 1 9 ^12 9 ^16 9 

plate plate plate 
*-*22 9 ^26 9 ^66 


, Plate f plate f plate 


D'„ , D 


D' 


12 9 ^ 16 


. plate . plate . plate 

r^' rv 


D 22 , D 


D' 


26 9 ^ 66 


dimension of hexagon- shaped stiffener pattern (see figures 21 
and 25), in. 

dimension of star-shaped stiffener pattern (see figures 23 
and 27), in. 

equivalent-plate coupling stiffnesses (see equation (22a)), lb 

plate contribution to equivalent-plate coupling stiffnesses 
(see equations (23b) and (44b)), lb 

contribution of sandwich face plate to equivalent-plate 
coupling stiffnesses (see equation (46b)), lb 

stiffener contribution to equivalent-plate coupling stiffnesses 
(see equation (23b)), lb 

dimension of hexagonal beam grid (see figure 21), in. 

beam-stiffener constitutive matrix defined by equations (A25) 
and (Dl) 

transformed beam constitutive matrix defined by equation (36) 
with elements given in Appendix D 

beam constitutive matrix defined by equation (32) 

stiffener spacing (see figures 4 and 5), in. 

equivalent-plate bending and twisting stiffnesses (see equation 
(22a)), in.-lb 

contribution of sandwich-plate core equivalent-plate bending 
stiffnesses (see equation (47c)), in.-lb 

plate contribution to equivalent-plate bending stiffnesses 
(see equation (23c) and (44c)), in.-lb 

contribution of sandwich face plate to equivalent-plate 
bending stiffnesses (see equation (46c)), in.-lb 
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stiffener 

D„ , 

stiffener 

D.2 , 

stiffener 

D,6 ■ 

stiffener 

stiffener 

stiffener 

C>22 , 

D’6 , 

d 66 


e 

e i’ e 2 


stiffener contribution to equivalent-plate bending stiffnesses 
(see equation (23c)), in. -lb 

sandwich face plate eccentricity (see figure 28), in. 

face plate eccentricities of sandwich plate (see figures 25-27), in. 


E, G 
E G 

J -'S’ '“'s 


Ex. E y , E z 


[E] 


beam membrane strain (see equations (11) and (A3)) 

extensional and shear moduli of isotropic materials, psi 

effective extensional and shear moduli of nonhomogeneous 
stiffeners (see equations (A13) and (A15)), psi 

principal extensional moduli for a homogeneous specially 
orthotropic material (see equation (A10)), psi 

strain-equivalence matrix defined by equations (30) 

beam- member strain energy (see equation (31)), lb-in 2 


S x strain energy density of equivalent plate (see equations (39)), lb 


/v P late 

S E contribution of plate wall to strain energy density of equivalent 

plate (see equations (39)), lb 

G xy , G yz , G xz principal shear moduli for a homogeneous specially 

orthotropic material (see equation (A10)), psi 


gL, g xz effective shear moduli for a nonhomogeneous beam 

(see equation (A13)), psi 

h plate-wall thickness of stiffened plate (see figures 5, 8-9, and 

figure 28), in. 

h c core thicknesses of sandwich plate (see figures 25-27), in. 

h s thickness of equivalent-stiffener layer (see figure 9), in. 

h , , h, face sheet thicknesses of sandwich plate (see figures 25-27), in. 


H 


dimension of star-shaped stiffener pattern (see figures 23 
and 27), in. 
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^YY’ ^ZZ’ ^ YZ 


moments and product of inertia of homogeneous beam 
(see equations (A8)), in 4 


lyY’ ^ZZ5 ^ > 


effective moments and product of inertia of nonhomogeneous 
stiffeners (see equations (A13)), in 4 


f 11’ ft’ I 




J 

Js 


effective moments and product of inertia of nonhomogeneous 
stiffeners with respect to the stiffness-weighted centroid (see 
equations (A28)), in 4 

torsional constant of homogeneous isotropic beam (see 
equation (A9)), in 4 

effective torsional constant for nonhomogeneous stiffeners 
(see equation (A15)), in 4 


k Y , k z 


inplane and transverse shear correction factors for homogeneous 
beams, respectively (see equations (A7)) 


k 


S 

Z 


j^plate 


L 

L s 

L , L 

x 7 y 

m y (x), m z (x) 


inplane and transverse effective shear correction factors for 
nonhomogeneous stiffeners (see equation (A24)), respectively 

transverse-shear correction factor for plates (see equation (44d)) 

hexagon side length shown in figure 21, in. 

generic beam-stiffener length (see equation (31)), in. 

beam-grid stiffener spacing (see figures 14-18,21, and 26), in. 

beam bending moments (see figure 11), in-lb 


7ft ^ %’ yy , equivalent-plate bending stress resultants in (x, y, z) global plate 

coordinates (see equations (18) and figure 6), in-lb/in. 

7ft xx , 7ft xx , 7ft xx equivalent-plate bending stress resultants in (X, Y, Z) beam 

coordinates (see figures 6 and 13), in-lb/in. 


plate plate plate 

^ xx ,3£ xy 


plate-wall bending stress resultants in (x, y, z) global plate 
coordinates (see equations (21)), in-lb/in. 




plate 
XX ? 




plate 
YY ? 




plate 

XY 


plate-wall bending stress resultants (see equations (l)-(9)) in 
beam (X, Y, Z) coordinates, in-lb/in. 
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stiffener stiffener stiffener 

«xx ,7ft YY Y 


xx 7 yy 7 xy 


^XX^YY^XY 


equivalent-stiffener-layer bending stress resultants 
(see equations (l)-(9)), in-lb/in. 

equivalent-plate membrane stress resultants in (x, y, z) global 
plate coordinates (see equations (18) and figure 6), lb/in. 

equivalent-plate membrane stress resultants in (X, Y, Z) beam 
coordinates (see figures 6 and 13), lb/in. 


^. plate ptate plate 

^ , 7t *7t 

*^xx ’ yy ’ *^xy 


plate-wall membrane stress resultants in (x, y, z) global plate 
coordinates (see equations (21)), lb/in. 


_ plate _ plate _ plate 

71 71 71 

XX ’ '^YY 5 't'XY 


plate-wall membrane stress resultants (see equations (l)-(9)) in 
(X, Y, Z) beam coordinates, lb/in. 


stiffener stiffener 
%XX , ®YV , 2? 


stiffener 

XY 


P(X) 

Q xz > Qyz 


equivalent-stiffener-layer membrane stress resultants 
(see equations (l)-(9)), lb/in. 

beam axial force (see figure 1 1), lb 

equivalent-plate transverse shearing stress resultants in (x, y, z) 
global plate coordinates (see equations (18) and figure 6), in-lb 


Qxz’ Qyz 


equivalent-plate transverse shearing stress resultants in 
(X, Y, Z) beam coordinates (see figures 6 and 13), in-lb 


plate plate 

Vxz ’ ''Cyz 


plate-wall transverse shearing stress resultants in (x, y, z) global 
plate coordinates (see equations (21)), lb/in. 


plate plate 

%'XZ ? nZyz 


plate-wall transverse shearing stress resultants in (X, Y, Z) beam 
coordinates (see equations (l)-(9)), lb/in. 


Q 


stiffener 
XZ ’ 


Q 


stiffener 

YZ 


equivalent-stiffener-layer transverse shearing stress resultants 
(see equations (l)-(9)), lb/in. 


Qu( z ) 

h fx t d , t s 
T(X) 


[T ], [TJ 1 


transformed stiffnesses of laminated-plate theory 
(see equations (44)), psi 

core wall thicknesses of sandwich plate (see figures 25-27), in. 
beam torque (see figure 11), in-lb 

strain transformation matrix and its inverse (see equations (17)) 


11 



[TJ, [TJ 1 stress transformation matrix and its inverse (see equations (16)) 

[T ], [T] 1 transverse shearing strain and stress transformation matrix and 

its inverse (see equations (16) and (17)) 

u(X), v(X), w(X) beam displacements of points on the reference axis (see 

equations (Al), in. 

U(X, Y, Z), V(X. Y, Z), displacements of beam material points (see equations 

W(X, Y,Z) (Al)), in. 


V Y (X), V Z (X) beam transverse shearing forces (see figure 1 1), lb 

(x, y, z) global plate coordinates (see figure 4), in. 


(x', y', l ) local coordinates of sandwich face plate (see figure 28), in. 

(X, Y, Z) noncentroidal beam coordinates (see figure 6), in. 


y s , y s 


Y, Z 



Yxz’ Yyz 
Y XZ’ Y YZ 


F 

1 XY’ 



stiffness-weighted eccentricities of nonhomogeneous stiffener 
(see figure 6 and equation (A13)), in. 

centroidal coordinates of a homogeneous beam with respect to 

the (X, Y, Z) beam coordinates (see equations (A8) and 
figure 6), in. 

stiffness-weighted eccentricities of nonhomogeneous stiffener 
(see figure 6 and equation (A13)), in. 

plate transverse shearing strains with respect to (x, y, z) 
coordinates (see equations (19) and figure 6) 

plate transverse shearing strains with respect to (X, Y, Z) 
coordinates (see equations (10) and figure 6) 

beam transverse shearing strains (see equations (11) and (A3)) 


s xx , £ yy , Yxy, Yxz, Yyz, £ zz plate strains with respect to (x, y, z) coordinates (see equations 

(17) and figure 6) 

s xx , £yy, Yxy, Yxz, Yyz, Ezz plate and beam strains with respect to (X, Y, Z) coordinates 

(see equations (10) and figure 6) 


£„, S yy , y xy plate membrane strains with respect to (x, y, z) coordinates 

(see equations (19) and figure 6) 
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o 

^xx’ 


Y> 


{e b } 

{e„} 


plate membrane strains with respect to (X, Y, Z) coordinates 
(see equations (10) and figure 6) 

vector of beam strains (see equations (30)) 

vector of plate strains referred to (x, y, z) global plate 
coordinates (see equations (34)) 




vector of plate strains referred to (X, Y, Z) beam coordinates 
(see equations (30)) 




eccentricities of nonhomogeneous stiffener (see equations 
(A27) and (A29)), in. 


4 S 9 h S 


eccentricities of nonhomogeneous stiffener (see equations 
(All) and (A29)), in. 


K 


O 

XX? 


o 



o 



plate bending strains with respect to (x, y, z) coordinates 
(see equations (19) and figure 6), in 1 


Kx X , Ky Y , Kxy plate bending strains with respect to (X, Y, Z) coordinates 

(see equations (10)), in 1 

cr xx , cr xy , a xz , a yy , a yz , a zz beam and plate stresses with respect to (x, y, z) coordinates 

(see figure 6), psi 


a.. 


T 



V V V 

V XY’ V XZ’ v YZ 




V V 

V ZX’ v z 


(£> o, O 


Cp x (X), Cp Y (X), Cp z (x) 


beam and plate stresses with respect to (X, Y, Z) coordinates 
(see figure 10), psi 

beam twisting strain (see equations (A2) and (11)), in 1 

inplane and transverse shear parameters (see equations (29)) 

major principal Poisson’s ratios for a homogeneous specially 
orthotropic material (see equation (A10)) 

minor principal Poisson’s ratios for a homogeneous specially 
orthotropic material (see equation (A10)) 

stiffness-weighted centroidal coordinates of nonhomogeneous 
beams and centroidal coordinates of homogeneous beams 
(see equations (All) and figures 6 and 11), in. 

dextral rotations of the beam cross-section about the X-, Y-, and 
Z-axes, respectively (see equations (Al)) 


13 



angle of diagonal stiffeners (see figures 14, 16-18, 21 and 23), 
degrees 


<F 

o o 

Xy> Xz beam bending strains (see equations (11) and (A2)), in 1 

V F S stiffener angle (see figures 4, 6, and 11), degrees 


Introduction 

Stiffened plates and shells are common structural forms used in the aerospace industry. For 
example, the Space Shuttle liquid-hydrogen tank utilizes a metallic "orthogrid" construction, as 
shown in figure 1. This type of structural arrangement has two families of uniformly spaced 
stiffeners that intersect at right angles. Stiffened metallic shells with this type of construction are 
often referred to as geometrically orthotropic because their overall equivalent extensional and 
bending stiffnesses are similar to that of a corresponding homogeneous orthotropic material when 
the stiffnesses of the shell skin and stiffeners are homogenized. A similar densely stiffened 
metallic orthogrid cylinder is shown in figure 2. Another common structural arrangement that is 
used for launch vehicles, known as an "isogrid," is shown in figure 3. This type of structural 
arrangement has three identical families of uniformly spaced stiffeners that intersect one another 
at a 60-degree angle to yield a geometrically isotropic shell wall. For this type of metallic 
stiffened shell, the homogenized extensional and bending stiffnesses are like that of an isotropic 
material. 

Approximating the overall behavior of stiffened plate and shell structures by using 
homogenized, "equivalent" stiffnesses and "effective" thicknesses has been in use for many years 
(e.g., see references 1-8). In recent times, these plate and shell stiffness representations still see 
utility in early stages of building-block analysis and design approaches that have evolved and are 
used extensively by industry. These simplified approaches are particularly useful for navigating 
the design space rapidly to identify the "first cut" of optimal preliminary designs, which are 
particularly important to the design of lightweight, high-performance launch vehicles such as the 
ARES V, previously under development by NASA. 

The earliest works that utilize equivalent stiffnesses to analyze the behavior of stiffened plates 
appear to be that of Huber. 1-3 Originally, this approach was driven by a lack of analysis methods 
that could simulate adequately the discrete nature of stiffened structures and the availability of 
analytical solutions for bending and buckling of homogeneous orthotropic plates. Similarly, the 
earliest corresponding works for stiffened shells appears to be that published by Flugge 4 in 1932 
and 15 years later by van der Neut. 9 As time passed and physical understanding matured, theories 
for analyzing stiffened structures by using equivalent stiffnesses and thickness continued to 
receive attention. For example, Dale and Smith 7 used equivalent stiffnesses to study the buckling 
behavior of compression-loaded sandwich plates with an orthogrid core. In 1946, Smith et. al. 8 
presented an improved formulation that accounts for variations in the neutral-surface position 
associated with local interactions between a plywood plate and a stiffener. Similarly, in 1947, 
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Pfluger 10 presented an improved theory and applied it to the buckling of stiffened plates. Pfluger’s 
work is considered an improvement to Huber’s much earlier work because it includes a more 
accurate treatment of the shearing stresses. One of the earliest works that uses an effective plate 
thickness for designing stiffened plates was given in 1948 by Gomza & Seide. 11 In their study, 
the effective plate thickness is obtained by adding the stiffener cross-sectional area, divided by 
the stiffener spacing, to the plate thickness. An early work that developed equivalent-plate 
stiffnesses for an isotropic corrugated sheet sandwiched between two isotropic flat skins was 
presented by Libove and Hubka 12 in 1951. Basically, formulas for equivalent-plate elastic 
constants were presented, for use with available symmetric-sandwich plate theories 13 ' 15 that 
include inplane dilatation and shearing, pure bending and twisting, and transverse shearing. 
Formulas for equivalent-plate elastic constants associated with coupling between inplane 
dilatation and pure bending and between inplane shearing and twisting were also presented in 
anticipation of the extension of sandwich plate theories to unsymmetric sandwich constructions. 
Experimental results were reported that were characterized as being in "close agreement" with the 
theoretical predictions for the bending and transverse-shear stiffnesses in the direction 
perpendicular to the corrugations and for the twisting stiffness. Later, in 1952, Benscoter and 
MacNeal 16 presented an equivalent-plate theory, based on first-order difference equations, for a 
thin multicell low aspect ratio, supersonic wing with straight spars and perpendicular ribs that 
includes transverse-shear deformations. Similarly, Horvay 17 presented an equivalent-plate 
formulation for a plate-like grid with a stiffener arrangement forming honeycomb cells in 1952. 

What may be the earliest work to provide expressions for equivalent-plate stiffnesses of 
isotropic plates with integral stiffeners was presented by Dow, Libove, and Hubka. 18 In this early 
work, 12 independent elastic constants were derived, following a laborious approach, that 
correspond to the A n , A p , A 22 , A 66 , B n , B I2 , B„, B 66 , D n , D,„ D 22 , and D 66 terms appearing in the 
contemporary constitutive equations for laminated-composite plates (see reference 19) that relate 
stress resultants to membrane strains and bending curvatures. The expressions given for the elastic 
constants were obtained by identifying the fundamental repeating element of the stiffened plate 
and then replacing each stiffener in the repeating element with a homogeneous orthotropic plate 
that is perfectly bonded to the skin of the stiffened plate. The strains in the repeating-element 
stiffeners are related to the corresponding plate strains and the strain energy of the repeating 
element is determined in terms of the equivalent-plate strains. Equations relating the equivalent- 
plate strains and the equivalent-plate stress resultants are obtained by differentiating the strain 
energy, which yield the equivalent-plate stiffnesses. Several expressions for an average or 
equivalent-plate thickness were also given that depend on the stiffener arrangement. A somewhat 
related study that focused on the torsional stiffness of orthogonally stiffened plates was presented 
by Crawford and Libove. 20 Around the same time period, Hoppmann and his colleagues 21 ' 24 
conducted experiments to determine the bending and twisting stiffnesses of orthogonally stiffened 
plates, and used these stiffnesses to perform calculations for static bending and vibration, based 
on an orthotropic -plate theory. In 1956, Huffington 25 published an analysis for determining the 
equivalent-plate stiffnesses for orthogonally stiffened plates without stiffener eccentricity, with 
respect to the plate skin. Later, in 1957, Bodner 26 analyzed buckling of ring- stiffened cylinders 
subjected to hydrostatic pressure by using an equivalent shell-wall thickness and an effective 
moment of inertia for the rings and shell wall combination. In his analysis, the rings were 
presumed to contribute only to the circumferential membrane and bending stiffnesses, and 
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buckling loads were obtained by treating the stiffened shell as an equivalent orthotropic shell that 
is nonhomogeneous in the radial direction. An effective length between rings was used in the 
calculation of the bending stiffness such that the maximum circumferential stress can be 
computed by using basic ring theory. The direction of ring eccentricity, with respect to the shell 
wall, is not included in the analysis. 

During the last fifty years, the equivalent-plate or shell approach was still being used as a first- 
approximation method to understand the behavior of eccentrically stiffened and sandwich plates 
and shells, and in simulating the behavior of massive beam-like, plate-like, and shell-like 
reticulated orbiting space structures whose discrete-member analysis was beyond the capability 
of the computing resources available at that time. 27 ' 138 For example, Stroud 30 presented a derivation 
of five elastic constants associated with pure bending and twisting of corrugation-stiffened panels 
in 1963. These constants were obtained by examining the force-deformation characteristics of a 
repeating corrugation-plate-skin cell. Results of experiments conducted to determine the elastic 
constants showed favorable agreement with the proposed theory. In 1964, Meyer & Bellefante 31 
presented equivalent-plate elastic constants for inplane dilatation and shearing and pure bending 
of a skin stiffened by an array of stiffeners that enclose equilateral triangles. These constants were 
also obtained by examining the force-deformation characteristics of a repeating cell. Equivalent- 
shell stiffnesses were presented by Sewall, Clary, and Leadbetter 32 in 1964 for cylinders with an 
orthogonal arrangement of rings and stringers sandwiched between inner and outer skins. In this 
work, the rings are presumed to make a negligible contribution to the longitudinal inplane and 
bending stiffnesses, the stiffness associated with coupling between longitudinal-circumferential 
dilatation and coupling between longitudinal-circumferential bending, the inplane shearing 
stiffness, and the twisting stiffness. Similarly, the stringers are presumed to make a negligible 
contribution to the circumferential inplane and bending stiffnesses. In contrast, the stringers are 
presumed to contribute to the coupling between longitudinal-circumferential dilatation, the 
coupling between longitudinal-circumferential bending, the inplane shearing stiffness, and the 
twisting stiffness. Block, Card, Mikulas, McElman, and Stein 37- 38 ' 40- 41 presented and applied 
equivalent- shell stiffnesses for orthotropic ring- stiffened corrugated cylinders and ring-and- 
stringer- stiffened cylinders in 1965 and 1966. In this series of papers, a relatively simple strain- 
energy approach was used to determine stiffnesses that include inplane dilatation and shearing, 
pure bending and twisting, and coupling between inplane dilatation and pure bending and between 
inplane shearing and twisting associated with stiffener eccentricity with respect to shell wall mid- 
surface. Transverse shearing stiffnesses were not considered. Also in 1966, Singer et.al. 42 
published stiffener expressions similar to those given by Block, Card, Mikulas, McElman, and 
Stein for cylinders with eccentric orthogonal rings and stringers. Their stiffness expressions are 
based on the presumption that the normal strains vary linearly in the skin and stiffeners and are 
equal at the point of contact between the skin and a stiffener. In addition, the stiffeners do not 
contribute to the inplane shearing stiffness, and the torsional stiffness is obtained by adding the 
torsional stiffness of the stiffeners to that of the skin. In 1968, Jones 49-51 presented stiffness 
expressions similar to those given in references 36-39, but for laminated-composite cylinders. 
This work appears to be the first to express the equivalent- shell stiffnesses in the terminology that 
has become standard practice for laminated composites in a large portion of the world. In 1969, 
Soong 52 ' 54 presented a derivation of the stiffness expressions similar to those given in references 
37, 38, 40, and 41, based on a strain-energy approach, but for orthotropic cylinders reinforced 
with a balanced pair of spiral stiffeners that make an arbitrary angle with the cylinder generators. 
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Attempts were made in the 1970s to refine the equivalent-plate stiffnesses for isotropic plates 
and shells with eccentric orthogonal stiffeners. For example, Cusens et.al. 58 presented an analysis 
that accounts for the stiffener contributions to the stiffnesses associated with coupling between 
inplane stretching and between anticlastic bending. In their analysis, the coupling effects are 
weighted by the size of the area of contact at the stiffener intersections. Similarly, Nishino et.al. 63 
modeled the interaction of the shearing stresses between a plate and an eccentrically located 
orthogonal stiffening grid and obtained equivalent-plate twisting stiffnesses that include the 
effects of inplane shearing deformations of the stiffeners. Likewise, Nemeth 81 derived equivalent- 
plate stiffnesses for single-layer grids made of beam members with inplane and out-of-plane 
(transverse) shear flexibility in 1979. Stiffness expressions were presented in this study for 
several grid configurations and in-depth comparisons with corresponding results obtained from 
discrete finite element models of the grids are given. In 1979 and 1980, Ko 82 ' 84 presented 
equivalent stiffnesses for an isotropic corrugated sheet sandwiched between two isotropic flat 
skins in a symmetric manner, and for a similar symmetric honeycomb-core sandwich plate. These 
stiffnesses are modified forms of the corresponding stiffnesses presented by Libove and Hubka 12 
in 1951, that account for corrugated walls with nonuniform thickness, and include the effects of 
transverse-shear deformations. This approach was extended by Ko to hat-stiffened panels in 
1991. 108 Later, in 1985, Reddy et.al. 55 presented equivalent stiffnesses for a circular cylindrical 
shell stiffened by an internal grid. In particular, symmetrically laminated shell walls stiffened 
with rings, stringers, and a pair of identical helical stiffeners making an arbitrary angle with the 
shell generators were considered. Expressions are given that include transverse-shear stiffnesses 
in addition to stiffnesses that account for inplane dilatation and shearing, pure bending and 
twisting, and coupling between inplane dilatation and pure bending and between inplane shearing 
and twisting associated with stiffener eccentricity with respect to shell wall mid-surface. Also in 
1985, Kolpakov 97 presented a method for determining the equivalent stiffnesses of elastic 
frameworks that includes an analysis of a planar beam gridwork with an overall negative-valued 
Poisson’s ratio. Deb and Booton 102 presented similar equivalent-plate stiffnesses in 1988 for shear 
deformable, orthogonally stiffened isotropic plates with eccentric stiffeners. Also, in 1988, Boot 
& Moore 103 presented a detailed list of the factors affecting the validity of equivalent-plate 
stiffnesses, and examined the importance of neglecting the contribution of the stiffeners to the 
coupling between biaxial stretching, for orthotropic plates. In 1989, Bunakov and Protasov 105 
presented equivalent-continuum stiffnesses for a pair of identical helical stiffeners with a 
rectangular-cross-section attached to a shell. These stiffnesses are based on a micropolar 
continuum model and include transverse-shear and bending stiffnesses of the beam members that 
are associated with deformations within the tangent plane at each point of the shell. In 1990, 
Won 106 presented equivalent-plate stiffnesses for isotropic plates reinforced with eccentric, 
regularly spaced pairs of oblique stiffeners. In his analysis, a uniform biaxial stress state is 
presumed to exist at the stiffener joints and is used to obtain a stiffener contribution to equivalent- 
plate stiffness associated with "Poisson coupling" between biaxial deformations. In addition, 
axial stresses in the stiffeners are presumed to develop that resist plate-like inplane shearing and 
twisting deformations, which leads to additional equivalent-plate stiffness contributions. Three 
coupled partial differential equations that govern the bending response are derived, in terms of the 
inplane and out-of-plane displacements, by minimizing the potential energy of the plate- stiffener 
system. An approximate "Huber-type" differential equation (see references 1-3) is also given, in 
terms of the out-of-plane displacement, for an orthotropic plate. This simplified equation, given 
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in terms of the out-of-plane displacement, represents bending about an "equivalent," inextensible 
neutral surface. The Huber-type equation is obtained by solving the equivalent-plate constitutive 
equations for the inplane stress resultants in terms of the membrane strains, and then substituting 
the resulting expressions into the equivalent-plate constitutive equations for the bending stress 
resultants. This step yields the bending stress resultants in terms of the out-of-plane displacement 
and second-order derivatives of the inplane stress resultants. The approximate "Huber-type" 
differential equation is obtained by substituting the equivalent-plate bending stress resultants into 
the remaining out-of-plane equilibrium equation and then neglecting the second-order derivatives 
of the inplane stress resultants. 

In 1993, Pshenichnov 112 published a monograph dealing with reticulated plates and shells, with 
an emphasis on single-layer plate-like and shell-like lattice structures in which the stiffeners are 
not eccentric with respect to the shell middle surface. The equivalent stiffnesses presented are 
based on a classical shell theory (no transverse shear flexibility) and are obtained by using tensor 
transformations to equate beam strains with corresponding shell strains and by equating shell 
stress resultants with transformed beam forces that are uniformly distributed across an equivalent 
shell wall. Although the analysis is based on a classical shell theory, an attempt is made to include 
the effects of stiffener bending in the tangent plane by expressing the beam shearing forces that 
act in the tangent plane in terms of the derivatives of the corresponding beam moments. These 
tangential beam moments are expressed in terms of the beam bending strain, in the usual way, but 
the beam bending strain is obtained in terms of the shell tangential displacements and strains by 
considering deformation associated with rotation about the unit vector normal to the middle 
surface. Although this approach captures tangential stiffener bending effects, the effects cannot 
be represented directly in terms of the shell strains and, as a result, do not enter into the equivalent 
stiffness expressions for plate-like and shell-like lattices. 

In 1995, Jaunky et. al. 116 117 presented a refined smeared-stiffener theory for grid- stiffened 
laminated-composite panels, based upon the earlier work presented by Smith et.al. 8 for a plywood 
plate with a single central stiffener that has a rectangular cross-section. The refinement presented 
in their work accounts for the variation of the neutral surface caused by interactions between the 
skin and the stiffeners. Results presented by these authors show more accurate equivalent-plate 
stiffness predictions for selected cases, but the analysis used to obtain the equivalent-plate 
stiffnesses is far more involved than the earlier, less refined approach used in earlier works such 
as references 37, 38, 40, and 41. Equivalent- shell stiffnesses for circular cylinders made of 
laminated-composite materials and with rings, stringers, and a pairs of identical geodesic 
stiffeners were presented by Gerhard et.al. 118 in 1996 (see pp. 56-75). These equivalent-shell 
stiffnesses were obtained by using strain transformation equations to express the stiffener strains 
in terms of the equivalent-plate strains, and by using force transformation equations to relate the 
force in each stiffener to the equivalent-shell stress resultants. Equivalent stiffnesses for 
laminated-composite flat plates and circular cylindrical shells stiffened by a grid of beams were 
presented by Chen and Tsai 119 in 1996. In their study, generally laminated walls stiffened with 
ribs, stringers, and a pair of identical diagonal stiffeners with an arbitrary orientation angle were 
considered. Grid- stiffness expressions are given that include out-of-plane (transverse) and 
inplane shear flexibility of the stiffeners and inplane stiffener bending in addition to the usual 
stiffnesses that account for inplane dilatation and shearing and pure bending and twisting. 
Similarly, Wodesenbet et. al. 133 presented an improved smeared-stiffener theory for isogrid- 
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stiffened laminated-composite cylinders in 2003. In this theory, the stiffness contributions of the 
stiffeners are obtained by relating the beam strains to the plate strains with the standard strain- 
transformation equations, and by expressing the shell stress resultants in terms of the beam forces, 
and their moments, that act on each stiffener component within a unit cell (repeating element). 
The equivalent-shell stiffnesses are then obtained by using a superposition of the stress resultant 
for the skin and the stress resultants for the stiffeners within the unit cell that is weighted by the 
corresponding volume fractions for the unit cell. 

The literature examined during the present study suggests that treating stiffened plates as an 
equivalent-continuum plate remains a useful design practice, provided the limitations of the 
theory are kept in mind. The surveyed results also indicate that, to a large extent, the derivations 
of equivalent-continuum plate stiffnesses have been ad hoc, the terminology is antiquated, and 
transverse-shear deformations are often neglected. Moreover, the criteria for defining an effective 
thickness of the equivalent-continuum plate appears to be unclear. The major objective of the 
present study is to present two systematic methods for deriving first-approximation, equivalent- 
continuum stiffnesses for eccentrically stiffened plates, based on a first-order transverse-shear 
deformation plate theory (e. g., see reference 139). The first method uses equilibrium and 
compatibility in a direct manner for plates reinforced by one or more families of rectilinear 
stiffeners. The second method is closely related and is based on using a basic, repetitive cell of 
the stiffened plate and then defining an equivalence between the strain energy of the basic cell and 
the corresponding equivalent plate. This method is particularly useful for stiffener arrangements 
that are not rectilinear. A second objective is to present a systematic set of equations that can be 
used to determine the equivalent thickness of stiffened plates. Toward these objectives, the 
analysis approach and details for the direct equilibrium-compatibility method are presented first 
for plates stiffened by one or more families of continuous rectilinear stiffeners. This first section 
includes discussions of how statical and kinematical equivalence is achieved between the 
stiffened plate and the equivalent-continuum plate, and equivalent-plate stiffnesses are given for 
several stiffening arrangements. Comparisons of the equivalent-plate stiffnesses with 
corresponding previously published results are also given. Next, the details of the energy- 
equivalence approach are presented, and the method is applied to obtain equivalent-plate 
stiffnesses for several stiffening arrangements. Additional comparisons of the equivalent-plate 
stiffnesses with corresponding previously published results and with the stiffnesses obtained 
herein by using the direct equilibrium-compatibility method are also given. Then, additional 
analysis is presented that shows how to obtain equivalent-plate stiffnesses for sandwich plates 
with two nonidentical, generally laminated face plates and a core made from several arrangements 
of beam stiffeners. Equations are also presented that show several different criteria that can be 
used to select an equivalent thickness of the equivalent-continuum plate, which may be needed to 
perform collateral design calculations. 


Direct Equilibrium-Compatibility Method for Rectilinear Stiffener Families 

Consider a perfectly flat plate of arbitrary shape, as depicted in figure 4a. Material points of 
the plate are located by the Cartesian coordinates (x, y, z), and z = 0 corresponds to the plate 
midplane. The plate is stiffened by a unidirectional family of identical prismatic beams that 
generally have nonhomogeneous cross-sections and that are equally spaced, as depicted in figure 
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4b. An example of the plate cross-section A-A indicated in figure 4b is shown in figure 5. In this 
figure, the stiffener spacing is denoted by d s and the effective nonhomogeneous-stiffener 

eccentricity, defined in Appendix A, is given by z s . The stiffeners are also presumed to be 
perfectly bonded to a laminated-composite skin with thickness h, and make an angle V I J S with the 
x-axis of the plate. Material points of the beam are located by the local, noncentroidal Cartesian 
coordinates (X, Y, Z), and Z = 0 also corresponds to the plate midplane. The geometric 
relationship between the plate and local beam coordinates is shown in figure 6. As indicated in 

Appendix A, the coordinates of the centroidal axis are given by (x, Y, z) , with respect to the 
(X, Y, Z) coordinate frame. Similarly, the coordinates of the stiffness-weighted centroidal axis of 
nonhomogeneous beams are given by (x, y s , z s ) , with respect to the (X, Y, Z) coordinate frame. 

The additional (§, q, 'C) coordinate axes shown in figure 6 are used to locate points of the beam 
with respect to the stiffness-weighted centroid. 

The stiffeners are modeled in the present study with the Timoshenko-type first-order shear- 
deformation beam theory presented in Appendix A, based on the local noncentroidal Cartesian 
coordinates (X, Y, Z). In addition, the beam material is presumed to be specially orthotropic with 
respect to the (X, Y, Z) coordinate frame, but may be nonhomogeneous in the cross-sectional 
planes. In particular, the principal axes of orthotropy are aligned with the (X, Y, Z) coordinate 
axes. Moreover, the effects of cross-sectional warping restraint associated with torsion of 
noncircular cross-sections are neglected in the kinematic equations and each cross section is 
presumed to warp in an identical manner. Beam constitutive equations are also presented in 
Appendix A for the nonhomogeneous specially orthotropic material in terms of effective 
engineering constants. Using a nonhomogeneous beam theory permits the modeling of tailored 
beam stiffeners such as those depicted in figure 7. Based on Appendix A, the nonhomogeneous 
beam stiffener is assigned an effective axial modulus E s , an effective shear modulus G s , cross- 

s s s 

sectional area A s , effective moments of inertia I YY and 1 ^, effective product of inertia I YZ , and 
effective torsional constant J s , in addition to the effective eccentricities y s , z s , y s and s s . 

Effective shear correction factors for the beam are denoted by k Y and k z . 

For the structural arrangement depicted by figures 4-6, the repetitive element shown in figure 
8 is used herein to represent the essential features of the structure. The stiffener in this repetitive 

element is positioned so that the effective eccentricity y s = 0 . Although a T-shaped stiffener is 

shown in figure 8, the stiffener may have a nonhomogeneous cross-section with an arbitrary 
shape. In the analysis that follows, the repetitive stiffened-plate element is modeled as the 
equivalent-plate element shown in figure 9, which consists of a laminated-composite wall bonded 
perfectly to a single, equivalent-stiffener layer. The reference surface of the equivalent plate is 
selected as the midplane of the plate wall, for convenience. The equivalent- stiffener layer is 
presumed to contribute axial stretching and bending stiffness, inplane shear stiffness, and twisting 
stiffness in the Y-Z plane, consistent with the classical continuum model of deformation. The 
approach that is followed in this section is to first establish an equivalence between the internal 
forces acting on the repetitive stiffened-panel element and the equivalent-plate element. This 
equivalence is referred to herein as statical equivalence and is based on the presumption that the 
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variation of stresses across the width of the equivalent-stiffener layer can be neglected for 
relatively small stiffener spacings. The beam stresses and corresponding resultant forces and 
moments, defined relative to the (X, Y, Z) coordinates, that are considered are shown in figures 
10 and 11, respectively. Similarly, the beam forces and moments acting on the repetitive 
stiffened-plate element and the corresponding distributed stress resultants acting on the 
equivalent-plate element are shown in figures 12 and 13, respectively. Then, an equivalence 
between the stiffener strains and the corresponding strains in the equivalent-stiffener layer is 
established, which is referred to herein as kinematical equivalence. The kinematical equivalence 
is based on the presumption that the strain at any point of a stiffener is identical to the 
corresponding strain at the corresponding point in the equivalent-stiffener layer of the equivalent- 
plate element. In addition, it is presumed that the variation of strains across the width of the 
equivalent-stiffener layer can be neglected for relatively small stiffener spacings. With these 
equivalences established, the stiffener-force contributions to the equivalent-plate stress resultants 
are expressed in terms of the equivalent-plate strains, and the corresponding equivalent-plate 
constitutive equations are determined in terms of the plate wall and stiffener properties. 

Statical Equivalence 

To establish statical equivalence between the repetitive stiffened-panel and equivalent-plate 
elements, consider the definition of axial stress resultant for the equivalent plate in the X-Y-Z 
coordinate system associated with the stiffener, given by 



( 1 ) 


where is the axial stress, h is the plate thickness, and h s is the thickness of the equivalent- 
stiffener layer shown in figure 9. The integration is partitioned to obtain 


where 
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Here, 


plate 

XX 


refers to the usual definition of the stress resultant for a plate that is given by the first 


integral in equation (2). Using the definition for the axial beam load, P(X), defined by equation 
(A5a) in Appendix A and applying it to the equivalent- stiffener layer of the equivalent-plate 
element, with the presumption that the variation of o xx across the width of the equivalent- 
stiffener layer is negligible, gives 


21 



p = 


I ( a xx dYdZ = d s | 

J J Ac J | 


T + h S 


a X v dZ 


Using equations (3) and (4) gives 


and equation (2) gives 
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where d s is the stiffener spacing (see figure 5). Following a similar process of partitioning the 
integrations and relating the two stiffener coordinate systems, the remaining stress resultants 
acting on the edge shown in figure 13 are expressed as 
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Likewise, applying the definitions of the beam forces V Y (X) and V Z (X) and beam moments M Y (X) 
and T(X) defined in Appendix A (see figure 1 1) to the equivalent-stiffener layer of the equivalent- 
plate element and neglecting the variation of the stresses across the width of the equivalent- 
stiffener layer gives 
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With these results, equations (7a) - (7d) yield 
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It is worth noting that the beam moment M Z (X), given by equation (A5f), vanishes when the 
variation of a xx across the width of the equivalent- stiffener layer is neglected, consistent with a 
plate theory based on a classical continuum model. 


Kinematical Equivalence 

The next step in the analysis is to establish kinematical equivalence so that the stiffener strains 
can be related to the appropriate plate strains in the equivalent-stiffener layer of the equivalent- 
plate element. The general expressions for the strains in a plate, modeled with first-order 
transverse-shear deformation plate theory, are given by (e. g., see reference 139) 
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These strains are defined with respect to the reference surface (plate wall midplane) of the 
equivalent-plate element shown in figure 13. Moreover, £ xx and £ YY are the extensional strains 

of the plate midplane and y XY is the corresponding inplane shearing strain. Likewise, k xx and 
k yy are bending strains of the plate midplane and k xy is the corresponding twisting strain. The 

symbols y xz and y YZ denote the transverse shearing strains. The kinematical equivalence is 
obtained based on the presumption that the strain at any point of a stiffener is identical to the strain 
at the corresponding point in the equivalent-stiffener layer of the equivalent-plate element. In 
addition, it is presumed that bending of the stiffener in the plane parallel to the plate midplane is 
negligible and, as a result, the variation of strains across the width of the equivalent-stiffener layer 
can be neglected. Furthermore, it is presumed that the eccentric stiffener contributes only half of 
the inplane shearing strain and half of the change in surface twist of the equivalent- stiffener layer. 
This presumption is rationalized by noting that the inplane shearing strain and the change in 
surface twist of a plate are deformation measures that have contributions associated with both 
cross-sectional faces of a differential plate element, and that the shearing and twisting 
deformations of the stiffener only act on one face of the equivalent- stiffener layer of the 
equivalent-plate element. Therefore, this last presumption represents an averaging of the stiffener 
contribution to the corresponding overall plate strains. Thus, the kinematical equivalence yields 
the following expressions 
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In equations (1 1), t is the change in twist of the beam, associated with torsion, and % Y and Xz 
are the beam bending strains associated with the changes in curvature in the X-Z and X-Y planes, 

o o 

respectively. The symbols T XY and r xz are the transverse shearing strains in the X-Y and X-Z 
planes, respectively. Substituting these expressions for the strains into the stiffener constitutive 
equations, given by equations (A19)-(A21) and (A25) in Appendix A, gives 

P = E s A s (e xx + Zs K xx) (12a) 
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Substituting equations (12) into equations (5) and (8) yields 
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Equations (13) give the contributions of the stiffeners to the pointwise equivalent-plate stress 
resultants, in the X-Y-Z coordinate system, in terms of the corresponding plate strains. However, 
another issue must be considered. In particular, the shearing stress resultant, ^ xy , and twisting 
stress resultant, of classical and first-order transverse- shear deformation plate theories are 
based on the symmetries = Tt and W - 7H, of the corresponding stress resultants acting 
on the cross-sectional faces of a differential plate element. The stiffener is presumed to provide 
no inplane shear stiffness and twisting stiffness to the equivalent-stiffener layer in the inplane 
direction parallel to its axis (e.g., see reference 30, p. 30). Therefore, with regard to the 
constitutive equations, equations (7a) and (7c) must be modified to yield 
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(15c) 


where it is noted that matrices associated with coupling between inplane and out-of-plane defor- 
mations are identical in both equations. It is important to also note that the stiffnesses occupying 
the third row and third column of the matrices in equation (15a) are associated with transverse- 
shear deformations of the beam member in a plane parallel to the plate midplane. As such, these 
stiffnesses are referred to herein as "inplane transverse-shearing stiffnesses." In contrast, the stiff- 
ness occupying the second row and second column of the matrix in equation (15c) is associated 
with transverse-shear deformations of the beam member in the direction perpendicular to the plate 
midplane and is referred to herein as an "out-of-plane transverse-shearing stiffness." 

The next step in the analysis is to relate the plate stress resultants in the X-Y-Z coordinate 
system to the x-y-z coordinate system shown in figure 6. Using the standard stress and strain 
transformation equations for the rotation of coordinates depicted in figure 6 gives 
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Based on the standard definitions of plate stress results, it follows that 
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From equations (10) and the corresponding equations in the x-y-z coordinate system, it follows 
from equations (17a) and (17b) that 
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with equations (18) it follows that 
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Substituting equations (19) into equations (15) and then substituting the results into equations (21) 
yields 


where 


A 

11 

A 12 

A 

16 

B 

11 

B , 2 

B ,6 

A 

12 

A 22 

A 

26 

B 

12 

B 22 

B 26 

A 

16 

A 26 

A 

66 

B 

16 

B 26 

B 66 

B 

11 

B ,2 

B 

16 

D 

11 

-D 12 

D ,6 

B 

12 

B 22 

b 

26 

D 

12 

D 22 

d 26 

B 

16 

b 26 

b 

66 

D 

16 

d 26 

D 6 6. 


I Qyz\ _ A 44 A 45 ( Yy z \ 

\Q„/'[a. 5 a !5 J v " 


A ii A 12 A , 6 
A 12 A 22 A 26 
A 16 A 26 A 66 


plate plate plate 

A a A 12 A 16 

plate plate plate 

^ 12 22 ^ 26 

plate plate plate 

A 16 A 26 A 66 


stiffener stiffener stiffener 

An A 12 A 16 

stiffener stiffener stiffener 

+ A 12 A 22 A 26 

. stiffener stiffener stiffener 


B,i B 12 B | 6 

B 12 B 22 B 26 
B i6 B 26 B 66 


plate plate plate 

^11 ^12 ^16 

plate plate plate 

^12 ^22 ^26 

plate plate plate 

^16 ^26 ^66 


stiffener stiffener stiffener 

B ii B , 2 B , 6 

stiffener stiffener stiffener 

+ b ]2 b 22 b 26 

_ stiffener stiffener stiffener 


D,, D, 2 D,6 
D|2 D 22 D 2 6 
D, 6 D 26 D66 


plate plate plate 

D D D 

J-^ll VJ 12 lj 16 

plate plate plate 

D D D 

^ 12 ^22 ^26 

plate plate plate 

D D D 

^ 16 ■ L/ 26 J - y 66 


stiffener stiffener stiffener 

Dn D,2 Di6 

stiffener stiffener stiffener 

+ d,2 d 22 d 26 

stiffener stiffener stiffener 

D,6 D,6 D 4 6 


29 



with 


A 44 A 4 

A, 


. stiffener 

stiffener 

. stiffener 

A„ 

A 12 

a 16 

stiffener 

stiffener 

stiffener 

A 12 

A 22 

A 26 

. stiffener 

stiffener 

. stiffener 

a 16 

A 26 

a 66 


plate plate 

A 44 -A- 45 

plate plate 

A,, A« 


[T 1 

A a 


+ 


stiffener stiffener 

-A 44 A. 45 

stiffener stiffener 

A A« 


E S A S 


0 

0 


0 

0 

0 


4cL 


(23d) 


0 


0 


k Y G s A s 


[t1 

8 


(24a) 


„ stiffener stiffener stiffener 

B 1 , B 12 B 16 

stiffener stiffener stiffener 

^12 ^22 ^26 

„ stiffener stiffener stiffener 

b 16 b„ b, 


9 26 


9 66 


T, 


E s A s z s 


0 

0 


0 

0 

0 


[t,; 


(24b) 


k y G s A s z s 
4cL 


D„ 

d 12 

D 

stiffener 

stiffener 


D,2 

E>22 

D 

stiffener 

stiffener 


d 16 

E>26 

D 


16 

stiffener 

26 

stiffener 


= \t T 

o 


E I s 


0 

0 


0 

0 

0 


0 


0 


1 G S J S 
4 d. 


t«; 


(24c) 


stiffener stiffener 

A 44 A 45 

stiffener stiffener 

A 45 A 


L 55 


t. 


E zG s A s 


X 


(24d) 


The constitutive terms with the superscript "plate" are the usual stiffnesses of the Reissner- 
Mindlin-type shear-deformation plate theory that are found in reference 139. The explicit 
expressions for the stiffnesses associated with the stiffeners are given by 
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is an inplane-shear-deformation parameter associated with the shear flexibility of the stiffeners in 
a direction parallel to the plate midplane, based on a first-order transverse-shear deformation 
beam theory, and 


S 

T z = 



(29b) 


is a corresponding transverse-shear-defonnation parameter associated with the shear flexibility of 
the stiffeners in a direction perpendicular to the plate midplane. The effects of inplane- shear- 
deformation and transverse-shear-defonnation of the stiffeners can be neglected in equations 

(25)-(28) by setting x Y = 0 and r s z = 0 , respectively. 


Expressions like those given by equations (25) and (27) have been given by Heki and Saka 50-56 
for homogeneous isotropic stiffeners with negligible inplane-shear stiffnesses and appear to be in 
agreement for the most part. Some discrepancies exist for the term associated with the twisting 
stiffness. However, the expressions for the stiffnesses associated with membrane anisotropy, 
equations (25c) and (25e), are in complete agreement with those of Heki and Saka. Expressions 
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like those given by equations (25) and (27) have also been given by Slinchenko and Verijenko 129 
for homogeneous isotropic stiffeners with negligible inplane-shear and torsional stiffnesses and 
are in complete agreement. Furthermore, expressions like those given by equations (25)-(27) have 
been given by Won 106 for homogeneous isotropic beam members with rectangular cross-sections 
of equal depth and negligible transverse-shear stiffnesses. The expressions given by Won are for 
a pair of oblique stiffener families and include higher-order effects associated with the interaction 
of the plate wall with the stiffeners at the stiffener intersections, which are neglected in a first- 
approximation analysis. When these higher-order effects are neglected, and only a single stiffener 
family considered, the expressions given by Won are in complete agreement with the 
corresponding expressions given herein. Pshenichnov 112 and Sumec 107 have also presented 
expressions like those given by equations (25) and (27) for homogeneous isotropic members with 
negligible transverse-shear stiffnesses and are in complete agreement. 


Analysis Applications 

The equivalent-plate stiffnesses given by equations (25)-(28) represent the homogenization of 
a single family of unidirectional, equally spaced, identical stiffeners that are oriented at an angle 
'F s with respect to the plate x-axis, as shown in figure 6. Inspection of the analysis reveals that 
the analysis is easily extend to plates stiffened with multiple families of stiffeners by simply 
adding the stiffness contributions of each family. Several examples of applying the analysis are 
presented subsequently for plates with the stiffener arrangements depicted in figures 14-18. In 
these figures, stiffeners aligned with the global plate x- and y-axes are referred to as stringers and 
ribs, respectively. Other stiffeners are referred to as diagonal braces or, simply, as diagonals. All 
stiffeners are eccentric with respect to the plate midplane, unless stated otherwise. The notation 
used for the material and section properties and orientation angle of each stiffener family in these 
examples is given in Tables 1-3. In these tables, the stiffener extensional modulus, shear modulus, 
eccentricities, moment of inertia, torsion constant, and transverse-shear correction factors refer to 
the corresponding effective quantities defined in Appendix A for a nonhomogeneous, specially 
orthotropic beam. Furthermore, quantities associated with the stringers, ribs, and diagonals are 
identified by the lower case subscripts or superscripts "s," "r," and "d," respectively. Wherever 
two families of diagonals are involved, the quantities associated with the two families are 
identified by the subscripts or superscripts "dl," and "d2. " 

Orthogrid plate with doubly braced bays. Consider a laminated-composite plate that is 
eccentrically stiffened with the arrangement of stringers, ribs, and braces shown in figure 14a. The 
pockets, or bays, formed by the skin, stringers, and ribs are diagonally braced with two 
nonidentical families of stiffeners. The spacings of the ribs and stringers are denoted by L x and 
L , respectively. The stiffness expressions are obtained from equations (23) and (25)-(28) by 
applying equations (25)-(28) to each family of stiffeners with the attributes given in Table 1. The 
resulting equivalent-plate stiffnesses are given in Appendix B. For this general stiffener 
arrangement, the stiffnesses in Appendix B include a broad range of anisotropies characterized by 
the A 16 , A, 6 , A 45 , D 16 , D 26 , B n , B 12 , B 16 , B 22 , B 26 , and B 66 equivalent-plate constitutive terms appearing 
in equations (22). It is worth noting that by neglecting the contributions of the ribs in the stiffness 
expressions given in Appendix B, the resulting expressions will then correspond to the diamond- 
shaped stiffener pattern shown in figure 15. 
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Several corresponding results have been published for stiffened plates and plate-like lattices 
with homogeneous, isotropic stiffeners. In particular, the stiffnesses given in Appendix B reduce 
to those given by Karmakar 88 and by Gerhard et. al. 118 (see pp. 57-75) for the special case of 
identical braces and members with negligible transverse-shear and torsional stiffnesses. 
Similarly, for the special case of no braces (orthogrids) and members with negligible transverse 
shear stiffnesses, the stiffnesses given in Appendix B reduce to those given by Block, Card, and 
Mikulas. 38 Likewise, the membrane stiffnesses given in Appendix B can be obtained from the 
Heki-Saka 50-56 equations that are similar to those given by equations (25), for the special case in 
which the stiffeners have negligible transverse-shear stiffnesses in the plane parallel to the plate. 
Expressions have also been given by Chen and Tsai 119 that include orthogrids and grids with only 
diagonal braces as special cases. The equivalent membrane and out-of-plane bending and 
transverse-shearing stiffnesses given by Chen and Tsai are in complete agreement with the 
corresponding results given in Appendix B. However, inplane transverse-shearing stiffnesses are 
presented by Chen and Tsai that disagree with the corresponding stiffnesses given in Appendix 
B. This disagreement appears to be associated with the fact that the stiffnesses given by Chen and 
Tsai are derived by using a method that accounts for combined bending and shearing of the 
stiffeners in the plane parallel to the plate. This combined stiffener bending-shearing action is not 
included in the first-approximation analysis presented herein. Similarly, equivalent-plate 
stiffnesses are given in reference 8 1 for rectangular orthogrids with two identical diagonal braces 
per bay and homogeneous, isotropic stiffeners. The corresponding stiffnesses presented herein 
for this special case are only in agreement if the terms associated with shear deformation and 
torsion of the stiffeners are neglected. At least part of the discrepancies, particularly that 
associated with inplane shear deformation of the stiffeners, appears to be associated with the 1/2 
appearing in equations (lie) and (14a) being neglected in reference 8 1 . It is important to point out 
that neglecting the 1/2 appearing in equations (lie) and (14a) results in a lack of the symmetry in 
equations (15a) and (15b). 

Bunakov and Protasov 105 presented equivalent-continuum stiffnesses for the stiffener 
arrangement shown in figure 14a, but without stringers and ribs and with identical diagonals with 
a rectangular cross-section. These stiffnesses are based on a micropolar continuum model of 
deformation and include transverse-shear and bending stiffnesses of the beam members 
associated with deformations in a plane parallel to the plate midplane. When the inplane bending 
stiffnesses of the stiffeners and the micro-rotation normal to the plate midplane are neglected, the 
stiffnesses given in reference 105, which include inplane and out-of-plane transverse-shear 
stiffnesses of the beam members, are in complete agreement with the stiffnesses given in 
Appendix B. It is important to point out that the expressions for the stiffnesses A 66 and D 66 were 
obtained from the micropolar theory by noting that both of the asymmetric inplane shearing strain 
measures and the asymmetric twisting strain measures used in reference 105 reduce to one half 
the inplane shearing strain and one half the twisting strain, respectively, when the micro-rotations 
normal to the plate midplane are neglected. The stiffnesses A 66 and D 66 were then obtained by 
relating the average of the asymmetric stress resultants to the corresponding strain. 

Orthogrid plate with singly braced bays. Equations (23) and (25)-(28) were also applied to 
a laminated-composite plate that is stiffened with eccentric stringers, ribs, and nonidentical 
diagonal braces, as shown in figure 16a. For this case, alternating bays, enclosed by stringers and 
ribs, are braced with a single diagonal stiffener, and the spacing of the ribs and stringers is denoted 
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by L x and L , respectively. The stiffness expressions are obtained from equations (23) and (25)- 
(28) by applying equations (25)-(28) to each family of stiffeners with the attributes given in Table 
2. The resulting equivalent-plate stiffnesses are also given in Appendix B. For the special case of 
homogeneous, isotropic stiffeners, identical braces, and members with negligible inplane 
transverse-shear stiffnesses, the stiffnesses given in Appendix B reduce to those given by Reddy 
et. a I." and by Chen and Tsai. 119 However, the out-of-plane transverse-shear stiffnesses of Reddy 
et. al. do not agree with those in Appendix B or with those given by Chen and Tsai, and appear to 
be in error. Likewise, for the special case of square bays with identical braces and members with 
negligible inplane and out-of-plane transverse-shear and inplane bending stiffnesses, the 
stiffnesses given in Appendix B reduce to those given by Hefzy & Nayfeh. 90 ' 99 

Isosceles-triangle stiffener pattern. The analysis of this section was also applied to a plate 
stiffened with an eccentric grid that has pockets that are isosceles triangles, as shown in figure 
17a. The notation used for the material and section properties and orientation angle of each 
stiffener family is given in Table 3, and the horizontal spacing of the triangles and the spacing of 
the stringers are denoted by L x and L , respectively. For this stiffening grid, the two families of 
diagonal beams are nonidentical and the stiffeners intersect at the vertices of the triangular 
pockets. The resulting equivalent-plate stiffnesses obtained by applying equations (23) and (25)- 
(28) are given in Appendix C. Expressions are also given in Appendix C for the special case of 
pockets that are equilateral triangles and members with identical properties. This particular 
stiffener configuration is often referred to as an isogrid because of the isotropic nature of its 
equivalent-plate inplane and bending stiffnesses. Similarly, equivalent-plate stiffnesses are given 
in reference 81 for this special case with homogeneous, isotropic stiffeners that are in complete 
agreement with those presented herein, except for the stiffness contributions associated with 
inplane shear deformation and the twisting stiffness D 66 . The discrepancy also appears to be 
associated with the 1/2 appearing in equations (lie) and (14a) being neglected in reference 81. 
The stiffnesses for this special case with homogeneous, isotropic stiffeners that have negligible 
transverse-shear stiffnesses are in complete agreement with the corresponding stiffnesses given 
by Hefzy & Nayfeh. 99 

Kagome stiffener pattern. Equations (23) and (25)-(28) were also applied to a plate stiffened 
with a similar eccentric stiffening grid presented by Wodesenbet et.al. 133 and shown in figure 18a. 
This type of stiffening configuration is sometimes referred to as a Kagome pattern. 140141 The 
notation used for the material and section properties and orientation angle of each stiffener family 
is the same as that given in Table 3. For this stiffener arrangement, the two families of diagonal 
beams are nonidentical and the horizontal stiffeners do not intersect the diagonal stiffeners at the 
triangle vertices. The resulting equivalent-plate stiffnesses are identical to those given in 
Appendix C for the isosceles-triangle stiffener pattern. For the special case in which the beam 
members are made of a homogeneous, isotropic material, have rectangular cross-sections, and 
negligible inplane transverse-shear and torsional stiffnesses, the stiffnesses given in Appendix C 
reduce to those given by Vasiliev et. al. 131 In addition, for the case in which all stiffeners have the 
same properties and all stiffeners have negligible transverse-shear, bending, and torsional 
stiffnesses, the stiffnesses given in Appendix C reduce to those given by Heki 50 and by 
Wodesenbet et.al. 133 Moreover, the stiffnesses given in Appendix C reduce to those given by 
Wang et. al. 138 for the special case of negligible transverse-shear and torsional stiffnesses. 
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Basic-Cell Energy-Equivalence Method 

In this method, a basic repeating cell of the stiffened plate, such as the cell shown in figure 
14b, is identified. The criterion for selecting the basic cell used herein is that it be a simple unit 
that will generate a stiffened structure of a given species by translation over the plate midplane, 
without any overlapping. Each basic cell contains a finite number of stiffener members. Consider 
a single, arbitrary beam stiffener. Applying the kinematical equivalence described by equations 
(11), the beam strains are expressed in terms of the equivalent-plate strains, in the (X, Y, Z) 
coordinates shown in figure 6, by 
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In these equations, the uppercase subscript, P, denotes reference to the (X, Y, Z) coordinates. The 
matrix defined by equation (30d) is referred to herein as the strain-equivalence matrix. By using 
these equations, the strain energy of the beam, given by (A25) and (A26), is expressed as 
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and where L s is the beam-member length. Performing the matrix multiplication gives 
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The elements of this matrix are based on the notation used in Appendix A for the effective 
engineering constants of a nonhomogeneous, specially orthotropic beam. In particular, the 
nonhomogeneous beam stiffener is assigned an effective axial modulus E s , an effective shear 

S S 

modulus G s , cross-sectional area A s , effective moments of inertia I YY and Izz, effective product 

S 

of inertia I YZ , and effective torsional constant J s , in addition to the effective eccentricities y s , z s , 
y s and i s . Effective shear correction factors for the beam are denoted by k Y and k) . 

Next, the plate strains in the (X, Y, Z) beam coordinates are defined in terms of the global plate 
(x, y, z) coordinates shown in figure 6 as follows 
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and where the lowercase subscript, p, denotes reference to the (x, y, z) coordinates of the plate. 
The transformation sub-matrices appearing in equation (34a) are defined by equations (17). In 
addition, 
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The beam strain energy defined by equation (31) is expressed in terms of the global plate 
coordinate system by 
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The explicit expressions for the elements of this matrix are given in Appendix D. 

The basic presumption of this energy-based homogenization method is that the basic repetitive 
cell is small enough, compared to the overall structural dimensions, so that the strain-energy 
density of every beam member within it can be represented adequately by a constant value. Based 
upon this discretization concept, the integrand of equation (35) is a constant and, as a result, the 
beam-member strain energy is given by 


= (37) 

The beam-member strain energy is now completely characterized by the plate strains and beam 
properties. The dependence on the strains is expressed functionally by 
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The "average" two-dimensional strain-energy density £, for the basic equivalent-plate cell, 
with area A u , is given by 
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This expression represents an equivalent pointwise function for the contribution of all beam 
members within the basic, repetitive cell to the equivalent-plate strain energy density, plus the 
contribution of the plate wall. The strain energy density in this expression is given in terms of all 
stiffener properties, stiffener orientations, and wall properties, in addition to the global plate 
strains. The next step in this method is to relate equation (39a) to the general form for an 
anisotropic plate, based on the first-order transverse-shear deformation theory presented herein. 
In particular, the strain-energy density of an anisotropic plate is given by 
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(40) 


By using the anisotropic -plate constitutive equations 
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the strain energy density given by equation (40) becomes a quadratic function of the plate strains. 
Therefore, it follows that the symmetric constitutive matrices are given by 
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Inspection of equation (39) indicates that the contributions of the stiffeners in the basic cell to the 
constitutive equations are given by 
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Using equations (43), the constitutive equations for the equivalent plate are given by equations 
(22) and (23). 


Analysis Application 

Several examples of applying the basic-cell energy-equivalence method are presented 
subsequently for plates with the stiffener arrangements depicted in figures 14-24. In these figures, 
stiffeners aligned with the global plate x- and y-axes are referred to as stringers and ribs, 
respectively. Other stiffeners are referred to as diagonal braces or, simply, as diagonals. All 
stiffeners are eccentric with respect to the plate midplane, and each family of stiffeners has 
different properties unless stated otherwise. The notation used for the material and section 
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properties and orientation angle of each stiffener member of the basic cell used in these examples 
is given in Tables 4-9. In these tables, the member designation denotes numbers associated with 
the beginning and end of the members comprising a basic cell, as shown in figures 14b, 16b, 17b, 
18b, 21b, and 23b. Also, the stiffener extensional modulus, shear modulus, eccentricity, moment 
of inertia, torsion constant, and transverse-shear correction factors, refer to the corresponding 
effective quantities defined in Appendix A for a nonhomogeneous, specially orthotropic beam. 
Quantities associated with the stringers, ribs, and diagonals are identified by the subscripts or 
superscripts "s," "r," and "d," respectively. Wherever two families of diagonals are involved, the 
quantities associated with the two families are identified by the subscripts or superscripts "dl," 
and "d2." For the examples presented subsequently, the Mathematica® computer program 142 
presented in Appendix E was used to obtain the equivalent-plate stiffnesses. 

Orthogrid plate with doubly or singly braced bays. The basic-cell energy-equivalence 
method was applied to the orthogonally stiffened plates with diagonally braced bays shown in 
figures 14a and 16a. The basic cell for the stiffener grid in figure 14a is shown in figure 14b and 
consists of four beam members with the basic-cell beam-member attributes given in Table 4. The 
basic cell for the stiffener grid in figure 16a is shown in figure 16b and consists of six beam 
members with the basic-cell beam- member attributes given in Table 5. The ends of each beam 
member are indicated in figures 14b and 16b by the filled red circular symbols. For both stiffener 
grids, the rib and stringer spacings are denoted by L x and L y , respectively, and the length of the 

diagonals is given by L d = \/lC + L“ . The dashed rectangle shown in figures 14a and 14b 
indicates the perimeter of the basic repetitive cell and has the area A ccll = L x L v . Similarly, the 
shaded rectangle shown in figures 16a and 16b indicates the basic repetitive cell, with the area 
A cell = 4L.E . Applying the analysis described in this section to these two basic cells yields the 
same equivalent plate stiffnesses given in Appendix B that were obtained by using the Direct 
Equilibrium-Compatibility Method presented herein. 

Isosceles-triangle and Kagome stiffener patterns. The basic-cell energy-equivalence method 
was also applied to a plate stiffened with the stiffener arrangements shown in figures 17 and 18. 
The basic cell for the grid in figure 17a is shown in figure 17b and consists of three stiffeners with 
the properties given in Table 6. Likewise, the basic cell for the Kagome grid in figure 18a is 
shown in figure 18b and consists of four beam members with the properties given in Table 7. The 
ends of each beam member are also indicated in figures 17b and 18b by the filled red circular 
symbols. For both stiffener grids, the horizontal spacing of the diagonals and the stringer spacing 
are denoted by L x and L , respectively. The dashed hexagon shown in figures 17a and 17b and 
the dashed rectangle shown in figures 18a and 18b indicates the perimeter of each basic repetitive 
cell and has the area A cell = L x L and 2L L , respectively, and the length of the diagonals is given 

by L = Vl' + L y . The shaded region surrounding each cell perimeter indicates how the cell 

relates to the surrounding structure, for clarity. In addition, figures 19 and 20 indicate how the 
basic cells in figures 17b and 18b are translated to generate the stiffened-plate structure, 
respectively. Applying the analysis described in this section to the basic cells with the attributes 
given in Tables 6 and 7 also yields the same equivalent plate stiffnesses given in Appendix C that 
were obtained by using the Direct Equilibrium-Compatibility Method presented herein. 
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Hexagon-shaped stiffener pattern. Another stiffened plate considered in the present study is 
shown in figure 21, in which the plate wall is reinforced with a hexagonal grid of stiffeners from 
three families. The basic cell for this structure is shown in figure 21b and consists of five 
stiffeners with the basic-cell beam-member attributes given in Table 8. The ends of each member 
are indicated in the figure by the filled red circular symbols, and the length of the ribs is different 
from the length of the diagonals. The dashed hexagon shown in figures 21a and 21b indicates the 
perimeter of the basic repetitive cell and has the area A cell = 2a(b + c). The shaded region 
surrounding the cell perimeter indicates how the cell relates to the surrounding structure. Figure 
22 shows how the basic cell in figure 21b is translated to generate the stiffened-plate structure. 
Applying the analysis described in this section to the basic cell yields the equivalent plate 
stiffnesses given in Appendix F. Similar membrane stiffnesses have been given by Heki 50 for a 
hexagonal stiffener arrangement in which all beam members forming a hexagon have the same 
length and are made of a homogeneous, isotropic material. The stiffnesses presented in Appendix 
F do not agree with those in reference 50. For the special case of homogeneous, isotropic 
members with the same length and identical properties and negligible inplane and out-of-plane 
transverse-shear stiffnesses, the stiffnesses given in Appendix F are in complete agreement with 
those given by Hefzy and Nayfeh." Similarly, equivalent-plate stiffnesses are given in reference 
81 for the special case of a homogeneous, isotropic single-layer grid comprised of members with 
the same length and identical properties. The stiffnesses presented herein are in agreement with 
those in reference 8 1 for this special case, except for the stiffness contributions associated with 
inplane shear deformation. This discrepancy also appears to be associated with the 1/2 appearing 
in equations (lie) and (14a) being neglected in reference 81. 

Star-shaped stiffener pattern. The analysis was also applied to a plate stiffened with the star- 
shaped stiffener arrangement shown in figures 23. For this case, the star shapes are formed from 
two oppositely directed isosceles triangles with coincident centroids. The basic cell for this 
structure is shown in figure 23b and consists of twelve stiffeners from three families with the 
basic-cell beam-member attributes given in Table 9. The ends of each member are indicated in the 
figure by the filled red circular symbols. The dashed hexagon shown in figures 23a and 23b 
indicates the perimeter of the basic repetitive cell and has the area A cell = 4BH/3, where B and H 
are the base and height of either isosceles triangle forming the star shape. The shaded region 
surrounding the cell perimeter indicates how the cell relates to the surrounding structure. Figure 
24 shows how the basic cell in figure 23b is translated to generate the stiffened-plate structure. 
Applying the analysis described in this section to the basic cell yields the equivalent plate 
stiffnesses given in Appendix G. 


Equivalent-Plate Stiffnesses for Sandwich Plates 

Inspection of the analysis presented herein reveals that the stiffnesses given in Appendices B, 
C, F, and G can also be used to obtain the equivalent-plate stiffnesses for sandwich plates with 
two nonidentical anisotropic face plates, such as those shown in figures 25-27. This task is 
accomplished by using equations (23) based on two plate members, neither of which are located 
at z = 0. The stiffeners forming the core are presumed to be made of a homogeneous ortho tropic 
or isotropic material. For convenience, the equivalent-plate reference plane is taken as the 
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midplane of the core layer formed by the stiffeners such that all core anisotropies associated with 
coupling between membrane and bending actions vanish. The contribution of the face plates to 
the equivalent-plate stiffnesses are obtained by applying the integral definitions of the plate 
stiffnesses, in accordance with figure 28. In particular, the stiffnesses of a face plate are given, 
with respect to the global plate (x, y, z) coordinates shown in figure 28, by 
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where h is the face-plate thickness, e is the face-plate eccentricity, Q^z) are the transformed 
lamina stiffnesses of first-order shear-deformation laminated-plate theory (see reference 139), 
and K plate is a transverse-shear correction factor. For the core, transverse-shear correction factors 
are used for the beam members. Next, the local (x', y', z') coordinates shown in figure 28 are 
used to obtain 


. plate 

. plate 

. plate 


plate 


, plate 


, plate 

A„ 

A 

12 

a 16 


A n 

A' 

12 

A' 

16 

. plate 

. plate 

. plate 


plate 


, plate 


, plate 

A 12 

A 22 

A 26 

= 

A 12 

A' 

22 

A' 

26 

. plate 

. plate 

. plate 


plate 


, plate 


, plate 

A 16 

A 26 

A 66 


A' 16 

A' 

26 

A' 

66 


(45a) 


plate 

plate 

^ plate 



, plate 


, plate 


, plate 


T-. 1 P' a,e 

_ , Plate 

_ , Plate 

B„ 

K 

K 


A' 

11 

A' 

12 

A' 

16 


B n 

B 12 

B',6 

plate 

plate 

plate 



, plate 


, plate 


, plate 


plate 

plate 

plate 

b 12 

B22 

B 26 

= e 

A' 

12 

A' 

22 

A' 

26 

+ 

B 12 

B 22 

B 26 

plate 

plate 

plate 



, plate 


, plate 


, plate 


^ , Plate 

plate 

-r^ t P ,3te 

® 16 

B26 

B 66 


A' 

16 

A' 

26 

A' 

66 


B 16 

B 26 

B 66 


(45b) 


44 



plate 

plate 

plate 



, plate 


, plate 


, plate 


plate 

plate 

plate 



, plate 

plate 

plate 

D„ 

d 12 

d; 


A' 

11 

A' 

12 

A' 

16 


B n 

B l2 

B l6 


D' 

11 

D 12 

D 16 

plate 

plate 

plate 

2 


, plate 


, plate 


, plate 


plate 

plate 

plate 



, plate 

plate 

plate 

D [2 

d 22 

D 2 6 

= e 

A' 

12 

A' 

22 

A' 

26 

+ 2e 

B 12 

B 22 

B 26 

+ 

D' 

12 

D 22 

D 26 

plate 

plate 

plate 



, plate 


, plate 


, plate 


, Plate 

_ , Plate 

T-. / P 1 ^ 



, plate 

__ , Plate 

, Plate 

D 16 

D 26 

D 6 6 


A' 

16 

A' 

26 

A' 

66 


B' 16 

B 26 

B' 66 


D' 

16 

D 26 

D' 66 


. plate 

. plate 


plate 

A IP 13 ’ 6 

a 44 

a 45 


A 44 

A 45 

. plate 

. plate 

— 


plate 

[a 45 

a 55 J 


[A' 45 

A' 55 


(45c) 


(45d) 


where z = z' + e and 
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Equations (46) represent the plate stiffnesses of a given face plate that are calculated in the usual 
way (see references 19 and 139) in which z' = 0 represents the plate midplane. With these defi- 
nitions, the equivalent-sandwich-plate stiffnesses are given by 
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The analysis of the present study was applied to the sandwich plates with two nonidentical 
anisotropic face plates shown in figures 25-27. For these examples, the eccentricities of the face 

plates, plates no. 1 and 2, are given by e, = - 4(h , + h c ) and e 2 = y(h 2 + h c ) , with respect to the 


46 



midplane of the core, respectively, where h c is the depth of the core and h, and h, are the 
thicknesses of plates no. 1 and 2, respectively. The sandwich plate shown in figure 25 has a 
regular hexagon-cell core composed of two member types, each with a rectangular cross-section. 
One member type is aligned with the y-axis of the plate, as shown in figure 25, and is referred to 
herein as a rib. The other member type makes an angle with the x-axis and is referred to herein 
as a diagonal. The core has the same layout as the hexagonal stiffener grid shown in figure 21a 

and with the dimensions A = TT , i and c = i . The geometrical properties of the face 

plates and the core are also shown in figure 25, and the section properties of the core are given in 
Appendix H. The core contributions to the equivalent-plate stiffnesses in equations (47) are given 
in Appendix H and are obtained by specializing the general stiffnesses given in Appendix F for a 
hexagonal stiffener arrangement. For the special case in which the face plates are made of the 
same isotropic material, the core is made of a second isotropic material, the thicknesses of the core 
members are identical (t s = t d = t), the face plates have identical thicknesses (h, = h, = h), and shear 
deformation is neglected, the stiffnesses in Appendix H reduce to those given by Kalamkarov 110 
(see p. 205), with the exception of the contribution of the torsional stiffnesses of the core 
members. 

The sandwich plate shown in figure 26 has a regular "orthogrid" core with the same layout as 
the stiffening arrangement shown in figure 14a, but without the diagonal braces. The geometrical 
properties of the face plates and the core are shown figure 26, and the relationship of these 
properties to the orthogonal stiffeners is given in Table 4 and Appendix I. The core contributions 
to equations (47) are given in Appendix I and are obtained by specializing the general stiffnesses 
given in Appendix B. Expressions are also given in Appendix I for the special case in which the 
two face plates and the core are made of three different homogeneous orthotropic materials. 
When these expressions are further reduced for face plates with identical thicknesses (h, = h 2 ) and 
shear deformation is neglected, the stiffnesses reduce to those given by Kalamkarov 110 (see pp. 
200-201), with the exception of the torsional stiffnesses of the core members. 

A sandwich plate with a star-shaped, or star-cell, core is shown in figure 26 which has the same 
layout as the stiffening arrangement shown in figure 23a. The relationship of the geometrical 
properties to the core stiffeners is given in Table 9 and Appendix J. Specializing the general 
stiffnesses in Appendix G yields the core contribution to equations (47) given in Appendix J. 
Expressions are also given in Appendix J for the special case in which the star shapes forming the 
core are composed of equilateral triangles with coincident centroids. 


Equivalent-Plate Thickness 

The equivalent-plate thickness, h + h s , shown in figure 9, is often needed to perform auxiliary 
design computations. However, for the analysis presented herein, the results reveal that this 
thickness is not defined uniquely. For example, enforcing equivalence between the area of the 
stiffener and the area of the equivalent-stiffener layer gives the equivalent-stiffener layer 
thickness as 
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In contrast, the stiffness-weighted first moment of stiffener area about the y-axis, defined by 
equation (A13e), is given by 



E x Z dYdZ = A s z s 


(49) 


Applying this equation to the equivalent-stiffener layer is generally complicated for an arbitrary 
nonhomogeneous stiffener. For the simpler case of a homogeneous stiffener, equation (49) gives 
the equivalence statement 


h s (h s + h) 


2A s z s 

d s 


(50) 


which is solved to obtain h s . Likewise, the stiffness-weighted second moment of stiffener area 
about the y-axis, defined by equation (A13f), is given by 



E x Z 2 dYdZ = f YY 


(51) 


Applying this equation to the equivalent-stiffener layer, for the simpler case of a homogeneous 
stiffener, gives 


h 3 + 3h h 2 + 3h; h = _3_i 

AA s ^ 2 s 4 11 s a Ayy 


(52) 


which is solved to obtain h s . Similarly, applying the formula from elasticity theory for torsion of 
a homogeneous, isotropic beam with a rectangular cross-section given in reference 144 (see p. 89) 
to the equivalent-stiffener layer gives the transcendental equation 


h 3 s d s 

3 


1 _ _64_ tanh | 


3jt d s 

v 2h s J 


= J S 


(53) 


Thus, equations (48), (50), (52), and (53) give four different values for the thickness of the 
equivalent-stiffener layer, but other criteria based on equations (A 13) exist. In particular, equation 
(48) results from equivalence based on membrane action, and equation (50) results from 
equivalence based on coupling between membrane and bending action. Equation (52) results 
from equivalence based on bending action, and equation (53) results from equivalence based on 
twisting action. The appropriate formula to use for h s should be based on the dominant physical 
characteristic response involved in a particular auxiliary computation. 
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Concluding Remarks 


A survey of studies conducted since 1914 on the use of equivalent-plate stiffnesses in modeling 
the overall, stiffness-critical response of stiffened plates and shells, and two detailed, 
comprehensive derivations of first-approximation equivalent-plate stiffnesses have been 
presented. The derivations are based on the Reissner-Mindlin-type, first-order transverse- shear 
deformation theory for anisotropic plates. The first derivation is based purely on static and 
kinematic equivalence between a stiffened plate and its homogenized equivalent. The second 
derivation is based on equivalence of the strain-energy density and is more amenable to 
complicated stiffener arrangements such as the hexagonal and star-shaped stiffening 
arrangements presented herein. In both derivations, the stiffener members are modelled as beams 
that are shear deformable within and transverse to the midplane of the plate, consistent with the 
classical continuum model of solid mechanics. In general, each stiffener is prismatic and may 
have a nonhomogeneous cross-section made from specially orthotropic materials. This general 
formulation allows the computation of equivalent-plate stiffnesses for stiffened panels such as 
those in which the stiffener caps are reinforced with high-strength, pultruded rods. Consistent 
with a first-approximation theory, inplane bending of the stiffeners and total compatibility 
between the plate skin and the stiffeners are neglected. Both methods presented herein have been 
shown to yield identical results. 

Equivalent-plate stiffness expressions, and a corresponding symbolic manipulation computer 
program, have also presented for several different stiffener configurations. These expressions are 
very general and exhibit the full range of anisotropies permitted by the Reissner-Mindlin-type, 
first-order transverse-shear deformation theory for anisotropic plates. The expressions presented 
in the present study were also compared with available, previously published results. For the most 
part, the previously published results are for special cases of the general expressions presented 
herein. These previously published results and are almost in complete agreement with the 
corresponding results contained herein and plausible reasons for the discrepancies have been 
given. Analysis has also been presented that extends the use of the equivalent-plate stiffness 
expressions to sandwich plates with nonidentical, anisotropic face plates. In addition, several 
criteria for obtaining an equivalent-plate thickness for stiffened plates have been presented. 
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Table 1. Stiffener attributes for an eccentric, orthogonal stiffening grid with two diagonal braces 
per bay (see figure 14) 


Stiffener 

family 

Spacing, 

d s 

Angle, >T S , 
degrees 

COStP s 

sin^ 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity, 

z s 

Stringers 

L , 

0 

1 

0 

A s 

Zs 

Zs 

Ribs 

L, 

90 

0 

1 

A r 

z r 

Z r 

Diagonals #1 

L X Ly 

L d 

<t> 

L x 

L d 

L y 

Lj 

A dl 

Zdl 

^dl 

Diagonals #2 

L X Ly 

L d 

00 

o 

e 

Ld 

Ly 

Lj 

A d2 

Zd2 

^d2 


Table 1. Concluded 


Stiffener 

family 

Effective 
moment of 

s 

inertia, I Y y 

Torsion 

constant, 

h 

Shear 

correction 

factor, k Y 

Shear 

correction 

factor, k| 

Effective 

modulus, 

E s 

Effective shear 
modulus, 

G s 

Stringers 

I s 

J s 

k s Y 

k) 

E, 

G s 

Ribs 

I r 

J r 

k r y 


E r 

G r 

Diagonals #1 

h, 

h, 

k y 

kz 

E dl 

G dl 

Diagonals #2 

c 

■C 

k y 

k? 

E,, 

G d2 






Table 2. Stiffener attributes for an eccentric, orthogonal stiffening grid with one diagonal brace 
per bay (see figure 16) 


Stiffener 

family 

Spacing, 

d s 

Angle, '\> s , 
degrees 

COS 1 ?;, 

sink's 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity. 

Stringers 

L 

0 

1 

0 

A 

Zs 


Ribs 

L 

90 

0 

1 

A 

Z r 

Z r 

Diagonals #1 

2L x L y 

L d 

<t> 

L x 

L d 

Ly 

L d 

A dl 

Zdl 

^dl 

Diagonals #2 

2L x L y 

L d 

00 

o 

e 

L x 

L d 

L , 

L d 

A d2 

Zd2 

Zd2 


Table 2. Concluded 


Stiffener 

family 

Effective 
moment of 

■ T S 

inertia, f YY 

Torsion 

constant, 

h 

Shear 

correction 

factor, k Y 

Shear 

correction 

factor, k) 

Effective 

modulus, 

E s 

Effective shear 
modulus, 

G s 

Stringers 

I s 

J 

k Y 

k) 

E s 

G s 

Ribs 

I 

J, 

k Y 

kz 

E r 

G r 

Diagonals #1 

h, 

h, 

k Y 

kz 

E dl 

G dl 

Diagonals #2 

c 


k Y 

kz 

E,, 

G d2 





Table 3. Stiffener attributes for an eccentric, isosceles-triangle stiffener pattern (see figure 17) 


Stiffener 

family 

Spacing, 

d s 

Angle, <P S , 
degrees 

COStP;. 

sin*P s 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity. 

Stringers 

C 

0 

l 

0 

A s 

z s 

! s 

Diagonals #1 

L X Ly 

L 

<t> 

Lx 

2L 

Ly 

L 

A al 

Zdl 

?di 

Diagonals #2 

L x Ly 

L 

00 

o 

e 

L x 

2L 

Ly 

L 

A d2 

z d2 

z d2 


Table 3. Concluded 


Stiffener 

family 

Effective 
moment of 

s 

inertia, I Y y 

Torsion 

constant, 

b 

Shear 

correction 

factor, k Y 

Shear 

correction 

factor, k| 

Effective 
modulus, E s 

Effective shear 
modulus, 

G s 

Stringers 

I s 

J s 

k s Y 

k) 

E, 

G s 

Diagonals #1 

b. 

b, 

k dl 

k y 

kz 

E dl 

G dl 

Diagonals #2 

c 

b. 

k d2 

k y 

k? 

E* 

G fi 
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Table 4. Basic-cell beam-member attributes for an eccentric, orthogonal stiffener pattern with 
two diagonal braces per bay (see figure 14b) 


Member 

designation 

Starting 
coordinates, 
(x, y) 

Ending 
coordinates, 
(x, y) 

Length, 

4 

Angle, W s , 
degrees 

COS 1 ?,, 

sin'F s 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity. 

1-2 

O 

Jb 

i 

L y , o 

L i 

0 

1 

0 

A 

Zs 


3-4 

S\™ 

i 

o' 

-TH 

o' 

4 

90 

0 

1 

A r 

z r 

Z r 

5-6 

1 

to j - ' 

1 

to j -1 

L x ^ ■ > 
2 ’ 2 

4 

<t> 

L s 
L „ 

Ly 

Ld 

A dl 

Zdl 

^dl 

7-8 

L x L y 

T’ 2 

Lx L y 
2 ’ 2 

4 

00 

o 

e 

_ 4 

Ld 

L, 

L d 

A d2 

Zd2 

^d2 


Table 4. Concluded 


Member 

designation 

Effective 
moment of 

s 

inertia, I Y y 

Torsion 

constant, 

J s 

Shear 

correction 

factor, ky 

Shear 

correction 

factor, k| 

Effective 

modulus, 

E s 

Effective shear 
modulus, 

G s 

1-2 

I s 

J s 

Ly 

k) 

E s 

G, 

3-4 

I 

J r 

k r y 

kz 

F 

G r 

5-6 

I dl 

J d , 

k dl 

k y 

kz 

4, 

G dl 

7-8 

4 

4 

k d2 

k y 

k d Z 

E d2 

G d2 






Table 5. Basic-cell beam-member attributes for an eccentric, orthogonal stiffener pattern with 
one diagonal brace per bay (see figure 16b) 


Member 

designation 

Starting 
coordinates, 
(x, y) 

Ending 
coordinates, 
(x, y) 

Length, 

4 

Angle, 

degrees 

cos^ 

sin'll 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity, 

4 

1-2 

-L,0 

L,0 

2L 

0 

1 

0 

A s 

Z S 


3-4 

0, -L 

’ y 

0, L 

y 

2L 

> 

90 

0 

1 

A 

Zr 

Z r 

5-6 

-L , -L 

x’ y 

L , L 

x’ y 

24 

<t> 

L x 

L d 

4 

4 

4, 

Zdl 

^dl 

5-7 

-L , -L 

x’ y 

L , -L 

x’ y 

2 4 

0 

1 

0 

A s 

Z s 


5-8 

-L , -L 

x’ y 

-L ,L 

x’ y 

2L 

> 

90 

0 

1 

A 

Z r 

f r 

7-8 

L , -L 

x’ y 

-L ,L 

x’ y 

24 

00 

o 

e 

4 

4 

4 

4 

A d2 

Zd2 

Zd2 


Table 5. Concluded 


Member 

designation 

Effective 
moment of 

• T S 

rnertia, f YY 

Torsion 

constant, 

T 

Shear 

correction 

factor, k Y 

Shear 

correction 

factor, k| 

Effective 

modulus, 

4 

Effective shear 
modulus, 

G s 

1-2 

I s 

J s 

k Y 

kz 

E 

G s 

3-4 

I 

J r 

k y 

kz 

L 

G r 

5-6 

h, 

4 

ky 

kz 

4i 

4, 

5-7 

I s 

j s 

ky 

kz 

E s 

G s 

5-8 

I 

j r 

ky 

kz 

F 

G r 

7-8 

4 

4 

k“ 2 

Ky 

kz 

4 

G d2 


V 1 


l; + l: 


and A , , = 4L L 

cell x y 
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Table 6. Basic-cell beam-member attributes for a plate stiffened with an eccentric, isosceles- 
triangle stiffener pattern (see figure 17b) 


Member 

designation 

Starting 
coordinates, 
(x, y) 

Ending 

coordinates, 

(x,y) 

Length, 

4 

Angle, 

degrees 

COS 1 ?,, 

sinTf, 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity, 

4 

1-2 

O 

Jb 

i 

^,0 

L i 

0 

1 

0 

A 

z s 

4 

3-4 

L x _Ly 
“X’ 2 

L x Ly 
4 ’ 2 

L 

<t> 

L x 

2L 

L y 

L 

A dl 

Zdl 

^dl 

5-6 

L s Ly 

L x Ly 

L 

e 

c 

oc 

L x 

Ly 

A., 




4 ’ 2 

4 ' 2 



2L 

L 


^d2 

^d2 


Table 6. Continued 


Member 

designation 

Effective 
moment of 

• T S 

inertia, 1 Y y 

Torsion 

constant, 

h 

Shear 

correction 

factor, ky 

Shear 

correction 

factor, k| 

Effective 

modulus, 

4 

Effective shear 
modulus, 

G s 

1-2 

I s 

J s 

k Y 

k) 

E, 

G, 

3-4 

4 

4 

k dl 

k y 

kz 

4, 

G,„ 

5-6 

4 

4 

k y 

k d z 2 

4 2 

G d2 


L= V L x + L y 
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Table 7. Basic-cell beam-member attributes for a plate stiffened with an eccentric Kagome 
stiffener pattern (see figure 18b) 


Member 

designation 

Starting 
coordinates, 
(x, y) 

Ending 

coordinates, 

(x,y) 

Length, 

4 

Angle, 

IF 

degrees 

COS 1 ?,, 

sin'Pj 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity, 

8 S 

1-2 

L x _4 
2 

L x L y 

T ’ 2 

L i 

0 

1 

0 

A 

z s 

8 S 

3-4 

L x L y 

2 ’ 2 

L X Ly 

2 ’ 2 

L 

0 

1 

0 

A s 

z s 

* s 

5-6 

L x T 

2 

L X T 

2 ’H 

L 

<t> 

L x 

2L 

Ly 

L 

A dl 

Zdl 

^dl 

7-8 

L x T 

2 ’ y 

L x r 

2 ’ y 

L 

00 

O 

e 

L x 

2L 

Ly 

L 

A d ; 

Zd2 

^d2 


Table 7. Continued 


Member 

designation 

Effective 
moment of 

s 

inertia, I Y y 

Torsion 

constant, 

J, 

Shear 

correction 

factor, ky 

Shear 

correction 

factor, k| 

Effective 

modulus, 

4 

Effective shear 
modulus, 

G s 

1-2 

I s 

J s 

k Y 

k) 

E s 

G s 

3-4 

I s 

J s 

k Y 

k) 

E, 

G s 

5-6 

I dl 

4 

k dl 

k y 

kz 

4, 

G dl 

7-8 

4 

4 

k d2 

k y 

k? 

E d2 

G d2 


L = V L x + L y 
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Table 8. Basic-cell-member attributes for a plate stiffened with an eccentric, hexagonal stiffener 
pattern (see figure 21b) 


Member 

designation 

Starting 
coordinates, 
(x, y) 

Ending 
coordinates, 
(x, y) 

Length, 

4 

Angle, 

W 

1 S’ 

degrees 

cosO s 

sinO s 

Area, 

A s 

Effective 

eccentricity, 

z s 

Effective 

eccentricity, 

4 


0 c 

a b + c 

L 


a 

b 




2-1 

’ 2 

2’ 2 

2 

-O 

L 

L 

A d2 

^d2 

^d2 


0 - c 

a b + c 

L 


a 

b 




2-3 

’ 2 

2’ 2 

2 

-180 + 0 

L 

L 

A dl 

Zdl 

^dl 


0 c 

o £ 



0 

1 




2-5 

’ 2 

’ 2 

c 

90 



A 

A 

Z r 


o £ 

a b + c 

L 


a 

b 




5-4 

■ 2 

2’ 2 

2 

<l> 

L 

L 

A dl 

^dl 

4 


0, $ 

a b + c 

L 


a 

b 




5-6 

2 

2' 2 

2 

180-0 

L 

L 

A d2 

^d2 

^d2 


Table 8. Continued 


Member 

designation 

Effective 
moment of 

s 

inertia, I Y y 

Torsion 

constant, 

4 

Shear 

correction 

factor, k Y 

Shear 

correction 

factor, k| 

Effective 

modulus, 

4 

Effective shear 
modulus, 

G s 

2-1 

I d2 

4 

k“ 2 

k Y 

k? 

4|2 

G d2 

2-3 

4 

4 

k Y 

kz 

4 , 

G dl 

2-5 

i 

j r 

k r Y 

kz 

L 

G r 

5-4 

4 

4 

ky 

kz 

4 , 

G dl 

5-6 

4 

4 

ky 

kz 

E d2 

G d2 


L = \/ a 2 + b" 
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Table 9. Basic-cell-member attributes for a plate stiffened with an eccentric, star-shaped stiffener 
pattern (see figure 23b) 





Table 9. Concluded 


Member 

designation 

Effective 
moment of 

• T S 

inertia, 1 YY 

Torsion 
constant, J s 

Shear 

correction 

factor, k Y 

Shear 

correction 

factor, k| 

Effective 

modulus, 

E s 

Effective 

shear 

modulus, 

G s 

1-12 

c 

Jdi 

k“ 2 

k ¥ 

4 2 

E d2 

G d2 

1-2 

c 

Id, 

kf 

k dl 

E d , 

G d , 

3-2 

I s 

J s 

k S Y 

kz 

E 

G s 

3-4 

c 

Jdi 

k Y 

kz 

E d2 

G d2 

5-4 

Id, 

Id, 

k Y 

kz 

Ed, 

G d , 

5-6 

I s 

J s 

k S Y 

kz 

E 

G s 

7-6 

Id! 

I. 

k“ 2 

k y 

kz 

E d2 

G d 2 

7-8 

Id, 

Id, 

k Y 

k dl 

E d , 

G d , 

9-8 

l 

J s 

k s Y 

kz 

E 

G s 

9-10 

Id! 

I. 

k d2 

k y 

k d z 2 

E d2 

G d 2 

11-10 

b, 

Id, 

k d ; 

k dl 

Ed, 

G d , 

11-12 

I s 

J s 

k s Y 

kz 

E 

G s 
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External tank 



Solid-rocket 

booster 


Orbiter 


(a) Space Shuttle 



(b) Internal structure of external ta nk 


Figure 1 . Space Shuttle and external ta nk structure. 
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Figure 2. Densely stiffened orthogrid cylinder. 
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Figure 3. Isogrid-stiffened cylinders. 
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(b) Family of beams within a "small" rectangular plate region 


Figure 4. Unidirectionally stiffened plate geometry. 
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Stiffener spacing, d s 



Figure 5. Cross-section A-A of a uniformly stiffened panel shown in figure 4b. 


Generic stiffener 


Stiffness-weighted 
centroidal axis 



y-axis 


Plate midplane 


x-axis 


X-axis 


Figure 6. Plate and stiffener coordinate systems. 
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Plate 

\ 



Figure 7. Examples of a nonhomogeneous stiffener cross-section. 



t 


d s 

Figure 8. Repetitive stiffened-panel element and coordinate systems. 


Plate 



Plate 

midplane 


Figure 9. Equivalent-plate wall. 
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Stiffness-weighted 
centroidal axis 


i Z-axis 



Figure 10. Beam stresses at an arbitrary cross-section. 
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z-axis 



Figure 11. Beam force and moment resultants at an arbitrary cross-section. 
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▼ < 


Z-axis 



Figure 12. Stiffener beam forces and moments acting on cross-section of repetitive 
stiffened-panel element. 



Figure 13. Plate stress resultants acting on cross-section of equivalent 
stiffened-panel wall. 
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Diagonal 

braces 



(a) Stiffener arrangement 



(b) Basic Cell 


Figure 14. Orthogonal stiffener pattern with two diagonal braces per bay and basic cell. 
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Figure 15. Diamond-shaped stiffener pattern. 


Diagonal 

braces 



(a) Stiffener arrangement 



(b) Basic Cell 


Figure 16. Orthogonal stiffener pattern with one diagonal brace per bay and basic cell. 
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Braces 



y-axis 



(b) Basic cell and adjacent structure 


Figure 17. Isosceles-triangle stiffener pattern and basic cell. 
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Braces 



y-axis 



(b) Basic cell and adjacent structure 


Figure 18. Kagome stiffener pattern and basic cell. 
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Figure 19. Generation of an isosceles-triangle stiffener pattern obtained by 
translating the basic cell defined in figure 17. 



Figure 20. Generation of a Kagome stiffener pattern obtained by translating the basic 
cell defined in figure 18. 
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(a) Stiffener arrangement 



x-axis 


Figure 21 . Hexagon- shaped stiffener pattern and basic cell. 
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Figure 22. Generation of a hexagonal stiffener pattern obtained by translating the basic 
cell defined in figure 21. 
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(b) Basic cell and adjacent structure 


Figure 23. Isosceles-star- shaped stiffener pattern and basic cell. 
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Figure 24. Generation of an isosceles-star-shaped stiffener pattern obtained by 
translating the basic cell defined in figure 23. 
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z-axis 

I L 


Plate no. 2 


y-axis 



(a) Cut-away view of sandwich plate construction 


z-axis 



(b) Plate cross-section 


Figure 25. Sandwich plate with nonidentical anisotropic face plates and a hexagon-cell core. 
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z-axis 

n 


Plate no. 2 


y-axis 



(a) Cut-away view of sandwich plate construction 


z-axis 



(b) Plate cross-section 


Figure 26. Sandwich plate with nonidentical anisotropic face plates and an orthogrid core. 
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(a) Cut-away view of sandwich plate construction 


z-axis 



(b) Plate cross-section 


Figure 27. Sandwich plate with nonidentical anisotropic face plates and a star-cell core. 
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z-axis 

A 


z'-axis 



Figure 28. Arbitrary face plate of a sandwich plate. 
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Appendix A 

Equations of First-Order Transverse-Shear Deformation Beam Theory 


Let (X, Y, Z) denote the coordinates of the material particles forming a beam stiffener, with 
respect to the noncentroidal coordinate system shown in figure 6. The displacement of any 
material particle in the beam is given by the X-, Y-, and Z-components U(X, Y, Z), V(X, Y, Z), and 
W(X, Y, Z), respectively. Consistent with Euler-Bernoulli and Timoshenko’s first-order 
transverse-shear-deformation beam theories, these displacement components are expressed as 


U(X, Y. Z) = u(X) - Ytp z (X) + Zcp Y (X) 

(Ala) 

V(X. Y. Z) = v(X) — Zcp x (X) 

(Alb) 

W(X, Y. Z) = w(X) + Yqp x (X) 

(Ale) 

where u(X), v(X), and w(X) are the displacements of the corresponding point on the beam stiffener 
reference axis, and cp x (X), cp Y (X), and cp z (X) are the dextral rotations of the beam cross-section 
about the X-, Y-, and Z-axes, respectively. For these kinematics, the cross-sectional planes are 
presumed to be planes of elastic symmetry; that is, the beam material is, at most, monoclinic with 
respect to the longitudinal X axis (see Lekhnitskii, 143 p. 271 and p. 293). The effects of cross- 
sectional warping restraint associated with torsion of noncircular cross-sections are neglected in 
these kinematic equations and each cross section is presumed to waip in an identical manner. 
Substituting these displacement expressions into the linear strain-displacement relations for a 
three-dimensional elastic solid gives 

8 xx (X, Y, Z) = e“ x (X) + Y Xz(X) + Z X>) 

(A2a) 

8 YY (X, Y, Z) = e^x, Y, Z) = y YZ (x. Y, Z) = 0 

(A2b) 

Y XY (x.Y,z) = r" Y (x)-Zx°(x) 

(A2c) 

Yxz(X.Y.Z) = r xz (X) + Yx"(X) 

(A2d) 

where 

e xx( X) = 

(A3 a) 

r x Y ( X ) = ^- ( Pz 

(A3b) 

r xz (X) = + Tx 

(A3c) 
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(A3d) 


T°(X) 


acp x 

ax 


X>) 


dqp Y 

dX 


(A3e) 



acpz 

ax 


(A3f) 


The symbol T (x) represents the change in twist of the beam, associated with torsion. The 

o o 

symbols X Y (X) and X'/A) represent the changes in curvature in the X-Z and X-Y planes, 

o o 

respectively. The symbols r xY (x) and F xz (X) represent the transverse shear deformations in the 
X-Y and X-Z planes, respectively. In first-order shear-deformation beam theory, the set u(X), 
v(X), w(X), cp x (X), cp Y (X), and cp z (X) are the primary unknowns. In classical Euler-Bernoulli beam 
theory, the transverse shear deformations are presumed negligible, which yields 


Ty= - 


3w 

ax 


<Pz = 


dv 

ax 


and eliminates two of the unknowns. In addition, equations (A3e) and (A3f) yield 


(A4a) 

(A4b) 


x;( x ) = 

Xz( x ) = ■ 


a~w 

ax 2 


2 

a v 
ax 2 


(A4c) 

(A4d) 


In the Euler-Bernoulli and first-order shear-deformation beam theories, all stresses are 
presumed to be negligible except a xx , a XY , and a xz . These stresses are shown in figure 10. The 
corresponding resultant internal forces are defined, with respect to the (X, Y, Z) coordinates, as 


P(X) = f f a xx dYdZ 

J J A s 

(A5a) 

V Y (X) = f f a XY dYdZ 

J J A s 

(A5b) 

V z (X)=ff a^dYdZ 
J J A s 

(A5c) 
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T(X) = |Ya xz - Zcr XY )dYdZ 


(A5d) 


M y (X)= ( 1 Z a xx dYdZ 

J J A s 

(A5e) 

M Z (X) = | f Y a xx dYdZ 

J J A s 

(A5f) 


where A s denotes the cross-sectional area of the beam stiffener. For these definitions, the 
corresponding positive valued moments are shown in figure 1 1 . 

The constitutive equations for a homogeneous, isotropic beam are obtained by substituting the 
beams strains and stresses into the constitutive equations for a general three-dimensional body. 
The result is 


o.xx = f^xx (A6a) 

(t xy — Gy XY (A6b) 

°xz = Gy xz (A6c) 

where E and G denote Young’s extensional modulus and the shear modulus, respectively. 
Substituting equations (A2) into (A6), and then substituting the results into equations (A5) gives 

P = EA s (e xx + Y)(x + ZX°) (A7a) 

V Y = k y GA s ( r° XY - z T°) (A7b) 

V z = k z GA s (r° xz + Y T°) (A7c) 

T = GA S k z Y rL-k Y Z r’ Y + GI P X° (A7d) 

M y = EA S Z e xx + EI yz X z + EI yy x y (A7e) 

M z — EA S Y e xx + EI ZZ x z + E I Y zX y (A7f) 
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where (Y, Z) are the coordinates of the cross-section centroid, k Y and k z are transverse-shear 
correction factors, and it is noted that 


( ( Y dYdZ = A S Y 

J J A s 

(A8a) 

j | ZdYdZ = A s Z 

J J A s 

(A8b) 

i | Z 2 dYdZ = I YY 
J J A s 

(A8c) 

| | Y 2 dYdZ = I zz 

J J A s 

(A8d) 

| | YZ dYdZ = I YZ 

J Ja s 

(A8e) 

JXj Yl+z! ) dYdz=i ' 

(A8f) 


The symbols I YY and I zz represent the area moments of inertia with respect to the noncentroidal 
(X, Y, Z) coordinate system, and I YZ denotes the corresponding product of inertia. The symbol I p 
denotes the polar moment of inertia, with respect to the origin of the coordinate system. To 
account for torsion of noncircular cross-sections, based on the St. Venant Theory of Torsion, 
equation (A7d) is expressed as 


T = GA, 


k z Y r - k Y Z R 


+ GJ T 


(A9) 


where J is the torsion constant of elasticity theory. 144 

A set of constitutive equations for nonhomogeneous beams, similar to equations (A7) and to 
those given by Allen and Haisler 145 (see pp. 171-187), are obtained as follows. First, the beam 
materials are presumed to be, at most, specially orthotropic with respect to the beam coordinate 
system. The material constitutive equations are given by 
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1 

V YX 

V ZX 

E x 

Ey 

~ E z 

V XY 

1 

V Z y 

E x 

E y 

~ E z 

Vxz 

V YZ 

1 

E x 

E y 

E z 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

J_ 

g y: 

0 

0 


0 

0 

0 

0 

j_ 

G v - 


0 

0 

0 

0 

0 

_L 

Gv. 



(A10) 


where E x , E Y , and E z are the principal Young’s moduli; v XY , v xz , and v YZ are the major Poisson’s 
ratios; v YX , v zx , and v ZY are the minor Poisson’s ratios; and G XY , G YZ , and G xz are the principal shear 
moduli. For a nonhomogeneous cross-section, these material parameters are functions of the Y 
and Z cross-sectional coordinates. Based on equations (A2) and the presumption that a YY , a zz , 
and a YZ are negligible, equation (A 10) reduces to 



i 

E 


A- 0 0 0 0 0 


0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 
1 


0 0 0 0 


Gv 


0 0 0 0 0 


0 

_L 

Gv, 



(Alla) 


or 


f \ 


o xx 


0 


/ o 

- 

\ o 

r 

°XZ 


\ °XY J 



E x 

0 

0 

0 

0 

0 


(e \ 

^xx 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 


0 1 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

Gxz 

0 


Yxz 

0 

0 

0 

0 

0 

Gxy 


y Y xy J 


(Allb) 


Substituting equations (A2) into (A1 1), and then substituting the results into equations (A5) gives 


P(x) = 


JI 


E x dYdZ e xx + 


| f E x Y dYdZ X ; + ff E x Z dYdZ 

A c A c 


(A12a) 


97 



V Y (x) = 


Si 


G xy dYdZ T 


--JI 


G xy Z dYdZ x 


V z (x) = 


( ( G xz dYdZ T xz + f f G xz YdYdZT° 

J J A s J Ja<. 


= Jl G “ YdYdZr - Jl G xy Z dYdZ r; v + 

JJ (G xz Y 2 + G XY Z : )dYdZT° 


M y (x) = | | E X Z dYdZ e xx + | f E x YZ dYdZ + f f E x z' dYdZ )(“ 

»/ »/ A c 1/ »/ A c 


M z (x) = 


Si 

Si 


E X Y dYdZ e^ x + f f E x Y^ dYdZ J&+ l ( E x YZ dYdZ £ 

A c c/ c/ A c 


By defining the effective quantities 


■if). 


E, = — I I E v dYdZ 


■ iff. 


G X7 II G xz dYdZ 


■iff. 


G XY = -E G XY dYdZ 


Si 


E x Y dYdZ 


A S E S 


Si 


E x Z dYdZ 


z c = 


A S E S 


I YY — 


Si 


E x T dYdZ 


(A 12 b) 
(A 12 c) 

(A 12 d) 

(A 12 e) 

(A 12 f) 

(A 13 a) 
(A 13 b) 
(A 13 c) 

(A 13 d) 

(A 13 e) 

(A 13 f) 
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(A13g) 


Izz — 


XX 


E x Y' dYdZ 


lyz — 


XX 


E X YZ dYdZ 


I P = 


XX 


Y 2 G xz + Z 2 G XY )dYdZ 


/T~s 

V GxzG, 


XX 


G XZ Y dYdZ 


AsGxz 


XX 


G xy Z dYdZ 


z, = 


A s Gv 


(A13h) 


(A13i) 


(A13j) 


(A13k) 


equations (A 12) are expressed as 


P — E s A s ^e xx + y s % z + z s % Y j 

(A 14a) 

V Y = k Y G XY A s ^r xY — z s X J 

(A 14b) 

V z = k z G xz A s ^r xz + y s X J 

(A14c) 

S — S — / S S ' S O 

G xz A s y s l xz — k Y G XY A s z s f XY + V G \/G xv I P X 

(A14d) 

M y = E s A s z s e xx + E s I yz X z + E s I yy )( y 

(A14e) 

M z = E s A s y s e xx + E S I zz^z EsIyz% y 

(A14f) 


S 

where E s is the effective extensional modulus of the nonhomogeneous beam stiffener; G xz and 
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s 

G xy are effective shear moduli; y s , z s , y s , and z s are stiffness-weighted first moments of area; 

S S S 

I Y Y and I zz are effective, stiffness-weighted moments of inertia; I YZ is the corresponding 

S 

stiffness-weighted product of inertia; and I P is a stiffness-weighted polar moment of inertia. It is 
noteworthy to indicate that for homogeneous orthotropic materials with the principal axes aligned 

with the beam axes E s = E x , G xz = G xz , G XY = G XY , y s = ys = y , z s = z s = Z 5 I YY = I YY , Izz = Izz , 

I yz = I Y z , and Ip = G XZ I ZZ + G XY I YY . To account for torsion of noncircular cross-sections based on 
the St. Venant Theory of Torsion, equation (A14d) is expressed as 

T = ^G^A^r^ - k Y G XY A s f s F XY + G s J s x (A15a) 

where J s is an effective torsion constant that based on elasticity theory, and generally depends on 
the material composition in addition to the cross-section geometry, and 

G s -V 5K7 (A 15b) 

For a stiffener with a rectangular cross-section of height and thickness r, made of a 

homogeneous orthotropic material, Leknitskii 143 gives an expression for the torsional stiffness that 
is equivalent to 


J s = 


G„Vt 3 4 


\G xy t 


%£| 

( G xz ^ 

f y 

D. 

1 

1 

Tt 5 o' 

iG XY j 

p= 1,2, 3, ... 

5 

P 


tanhl 


[ i 

( G V* 

'-’XY 

^ G xz j 


pjt<^ 


2 1 


(A16) 


The strain energy of a beam is obtained from the general definition given by the theory of 
elasticity; that is, 


-i/JI (tfxxExx + GyyEyy + a zz£zz + °xyYxy+ °xzYxz + °yz Y Y z) dXdYdZ (A17) 


where V s denotes the volume of the material particles forming the beam stiffener. For the 
prismatic beams considered herein, the volume is the product of the cross-sectional area A s and 
the length L s . Using equations (A2) and (A5) with (A17) gives the beam strain energy as 

s z = \f (peL + M zX ; + m yX ; + v Y r XY + v z rL + Tx°)dx (Ai8) 

Jl s 

For the analysis presented herein, it is convenient to define a vector of beam forces and moments 
by 
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and strains by 


(A19) 


{?} = { p v Y v z m 2 m y t } 

{e b } =| e xx r XY rL Xz x°y t ° } (A20) 


where the superscript "T" denotes matrix transposition. Then, equations (A7) and (A9) yield 


where 


[eY 


EAc 


{?} = [<?]{«.} 


k Y GA s 


0 

0 

k z GA s 


symmetric 


EA S Y 

0 

0 

EI 77 


(A21) 


EA S Z 

0 

0 

FT 

FT 

- L ' A YV 


- k Y GA s Z 
k z GA s Y 
0 
0 

GJ 


(A22) 


is the constitutive matrix for a homogeneous beam made of isotropic material. For a 
heterogeneous beam made of, at most, specially orthotropic materials, equations (A14) and (A15) 
yield 



E S A S 0 0 

• k Y G XY A s 0 

• • k z G xz A s 

• • • 

• • • 

symmetric • • 


s A s y s E s A s z s 0 


0 

0 

k Y G XY A s z : 

S = 

0 

0 

kzGxzAsyj, 

E I s 

J ^S A ZZ 

E I s 

A “'S A YZ 

0 

• 

E I s 

A “'S A YY 

0 

• 

• 

G S J S 


(A23) 


To obtain a matrix equation similar in form to equation (A22), it is convenient to introduce the 
stiffness-weight shear correction factors 



k 


S 

Z 


= k z 



(A24a) 

(A24b) 
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such that equation (A23) becomes 



E S A S 0 0 E s A s y s E s A s z s 0 

• kyG s A s 0 0 0 -k S y G s A s z s 

• • k Z G s A s 0 Ok Z G s A s y s 

• • • E S I ZZ E s I yz 0 

• • * * E S Iyy 0 

symmetric • • • • G S J S 


(A25) 


Note that for a homogeneous beam, equation (A25) has the same structure as equation (A22). By 
using equations (All) - (A15), the strain energy is expressed as 



(A26) 


In equation (A25), (x, y s , z s ) define the coordinates of a stiffness-weighted centroidal axis, 

shown in figures 6 and 11, which is coincident with the centroidal axis for prismatic beams with 
homogeneous cross sections. By defining the stiffness-weighted-centroidal coordinates q, 'Q 
shown in figures 10 and 11, the effective quantities defined by equations (A 13) yield 


E s = -j— f f E x dqd^ 

Ja s 

(A27a) 

I G xz dqd^ 
A sJ Ja s 

(A27b) 

G xy = ~r f f G *y dqdc 
A s J J A s 

(A27c) 

f f E x q dqdC 

J J A s 

q s = =0 

a s e s 

(A27d) 

f f E x dqd'C, 

J J A s 

= =0 
a s e s 

(A27e) 
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S S _ 2 

I YY = Ipr) 

(A27f) 

S S 2 

Izz = Iix + A s y s 

(A27g) 

Iyz = + A s y s z s 

(A27h) 

y s = ^s + ys 

(A27j) 

Z s — z s 

(A27k) 


with the stiffness-weighted-centroidal moments and product of inertia 



E x C 2 dpd'C 

~E S 

E x rf dpd'C 

~Es 




E x pC dpdC 


E s 


and the shear-stiffness-weighted first moments of area 



(A28a) 


(A28b) 


(A28c) 


(A29a) 


(A29b) 


For beams with homogeneous cross sections, the shear-stiffness-weighted first moments of area 
vanish and the stiffness-weighted-centroidal moments and product of inertia become the 
corresponding centroidal moments and product of inertia. 
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Appendix B 

Equivalent-Plate Stiffnesses for a Plate Reinforced with Eccentric Stringers, 

Ribs, and Diagonal Braces 

The expressions presented in this appendix are for a general laminated-composite plate that is 
stiffened with stringers, ribs, and diagonal braces, as shown in figure 14 or figure 16. The 
stringers and ribs are eccentric with respect to the plate midplane and the pockets formed by the 
skin, stringers, and ribs are diagonally braced with two nonidentical families of eccentric 
stiffeners. The notation used for the material and section properties and orientation angle of each 
stiffener family is given in Tables 1 and 2 for the stiffener arrangements shown in figures 14 and 
16, respectively. In particular, the subscripts and superscripts "s" and "r" refer to quantities 
associated with the stringers and ribs, respectively. The subscripts and superscripts "dl" and "d2" 
refer to quantities associated with the two nonidentical families of diagonals. The stiffness 
expressions are obtained from equations (23) and (25)-(28) by applying equations (25)-(28) to 
each family of stiffeners with the attributes given in Tables 1 and 2. The resulting equivalent-plate 
stiffnesses are given as follows where a - 1 and 1/2 for the stiffener arrangements shown in 
figures 14 and 16, respectively. In the expressions that follow, the stiffener extensional modulus, 
shear modulus, eccentricity, moment of inertia, torsion constant, and transverse- shear correction 
factors refer to the corresponding effective quantities defined in Appendix A for a 
nonhomogeneous, specially orthotropic beam. In addition, the length of the diagonal stiffeners is 

given by L d = »Jl\ + L l ■ 


plate E,A, <*L X 

A„ = A„ +-^ + ^- 

L y L d L 


E d iA dl | 

i ,n 

l+xt'tf 

+ E d2 A d2 | 

i ,n 

l+xftf 


L - 


l L J 


(Bl) 


plate X L y T 

A ,2 = A 12 H 3 


E dl A dl (l -ty) + E d2 A d2 (l -t“) 


(B2) 


. A P laK 

A 16 = A 16 




f dl 

T v 

Ll 


( ty 

L y 

A 

+ 3 

E d iA dl 


1-^ 

- F A 

J ^d2- rv d2 


1-^ 


L d 



[ L X J 


l 

[ L X J 

J 


(B3) 


_ Plate EA, « L y 

^-22 — -^22 ~ t ~ 3 

L * L d L x 



( T A 

i+x d ;h 


( 2 A 

F A 1 

J -'dl- rv dl 

| + E d2 A d2 j 

l+xftf 


l L J 


l m] 


(B4) 


. plate y 

A 26 — A 26 h — 



f 

dl 

r 

T 2 ' 


f 

d2 

r 

T 2 " 


E d iA dl 

1 

_ 

1 

_ L x 

-FA 

J ^d2- rv d2 

1 

_ 

1 

_ L x 

2 



l 

2 

- 

L y J 

/ 

{ 

2 

- 

LyJ 



(B5) 
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A — A P ' a,e _ 1 _ 

66 “ 66 4 L„ 4 L V 


<*E X L T dl L L 

1+ tl!-l: 


r d2 T L 

+ E d ,A d , 1 + — E — 

4 L y L x 


3 / 2 \ / 2 

plate E„A„Z_ #L X _ dl= E l I d2= L 

B 1 1 — B u H _ I 5 E dl A dl z dl + x Y z d , , + E d2 A d2 z d2 + x Y z d2 j - 

L > L d L v L I L x 


n p ia'e ^L x L y r I dl= \ „ . /■_ d2= 

^ 12 — ^ 12 3 — Z dl J + Ed2^d2^ Z d2 — Z d2 

Lh 


_ plate ^ x T 7 ., L v 

B I6 = B; + — E dl A d , | _ y 

Ed \ E x 


<12= T - 

fa 7 - TYZd2 1 - y 

I^d2 rt d2 Z; d2 ^ 1 2 

Z L, 


” _t3 22 


plate EjA^j «L y | _ dl= L I „ . I _ d2= 

77 “I” ^ H” ^ Ehi^hi Z hi + T y Z H | I + Zj 7 + t y z 


dl^dl ^dl ^ L Y ^dl 2 ■ - L ^d2- r ^-d2 ^d2 ^ U Y ^d2 2 

L y L y 


/ dl = / T 2 

p _ p plate I y F A ~ T Y Z dl 1 ^ 

^26 ^26 ■" 3 L^dl-^-dl Z dl ~ * 2 


d2= 2 

E A — - -1— ( y 

H 7 A ,17 Z .|7 _ 1 3 


plate E S A S X Y Z S E r A r T Y Z r 

b 66 = b 66 +^ 7 — + ^— + 


#L X L dl = l L T d2 f L L 

l_ FA 7 I X Y Z dl , F A 7 . X Y Z d 2 ^x _ 

3 -^dl^dl Z dl 'a t T ^d2^d2 Z d2 ‘a t t 

L d 4 L y LI 4 L y L x 


E I 


^ nlate Hi, r ^x \j 

Du = Du "I" . I 3 E dl I dl + E d2 I d2 + G dl J d , + G d2 J d2 

\ . T T /I 


'y L d L y 


plate ^L x L y -1 

D 12 = D 12 h 3 E dl I dl + E d2 I d2 — G dl J dl — G d2 J d2 
L d J 
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(B15) 


plate L. 

Di 6 — D ]6 + Vt 3 

Lrl 


Edildi E d2 I d2 + 


G d2 J d2 G dl J d i L y 

o 1 - 7T 


D = D Pkte + — + ^ L> 

^22 J - y 22 ' T ' 3 

L * L x L d 


Edlldl + ^d 2 ^d 2 + (^dlJdl + 


plate L y 

U 2 6 — E »26 + # ; 


„ „ „ „ GdJdi G d 2 J d2 1 L x 

P | F T -i_ 1 1 

^dl^dl - L 'd2 A d2 ■ n I _ 2 1 

J 


D, 


plate 


lloj, . GJ.1 . gL - L > 


= D “ + *|l“ + Tf 1| + 

\ ‘-‘y ^x 


(E d il d i + E d2 I d2 j + 


#L x L y 


G J + G J " ^ 


3 l °dl J dl * VJ d2 J d2 II T T 

4L;, 1 ; \Ly L 


plate E r A r X z ^E y 


A 44 = A r 4 7 + + — — (E dl A dl x z + E d2 A d2 xI 2 

E x E x Iv d ' 


plate [ E dl A dl X z E d 2 A d 2 X z ^ 

A 45 = a 45 +j 




J 


plate 


EAxl ^ L x 


A 55 — A 55 + I' . . I x^dl-n-dl “Z * iJ d 2 - r *' d2 *'Z 

-^y -^y-^d V 


■j—i * dl -j— i » d2 

E H 1 A rtl T 7 + E H 9 A H 9 t 7 


where the inplane- shear-deformation parameters are defined as 


s k Y G s 


Ty " e. 


r k Y G r 

Tv = — — 


di k Y G dl 


Xv = 


d2 _ k Y G d2 

V V — „ 


(B 16 ) 

(B17) 

(B18) 

(B19) 

(B20) 

(B21) 

(B22) 

(B23) 

(B24) 

(B25) 
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and the transverse-shear-deformation parameters are defined as 


k,G 

T 7 = ^ 

E s 

(B26) 

k,G. 

T z = — — 

E r 

(B27) 

dl 

di _ k z G dl 

z “ E 

^dl 

(B28) 

d2 

d2 _ k z G d2 

T z - 

J -'d2 

(B29) 


Setting these shear-deformation parameters equal to zero eliminates the effects of stiffener shear 
deformation. For the special case of no diagonal braces, the stiffnesses reduce to 


A 11 A 12 A 

A i , A ,2 A 
A 16 A 26 A 
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26 
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12 
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E ; As o o 
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o E / Ar 0 

Q Q E s A s X Y 

4L y 4 L x 


(B30) 


B n B ,2 B 16 
B 12 B 22 B 26 
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16 

26 
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o E : Ar z. 
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n E A t y= t y _ 

(J z. H z. 

4L, s 4L V r 


(B31) 


D„ D, 2 D, 6 

Dj 2 D 22 d, 6 
d, 6 d 26 d 66 
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D 
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D 
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D 
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D 
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D 

plate 

16 

26 
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0 0 


0 

0 

l/M Gjl 

4 U, L.J 


(B32) 
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A 
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44 

A 45 
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45 
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E r A r Tz q 
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E s Ajz 
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Appendix C 

Equivalent-Plate Stiffnesses for a Plate Reinforced with an Isosceles-Triangle 

or Kagome Grid Stiffener Arrangement 

The expressions presented in this appendix are for a general laminated-composite plate that is 
stiffened with stringers and two nonidentical families of diagonal braces, as shown in figure 17 or 
figure 18. The stringers and diagonal braces are eccentric with respect to the plate midplane and 
the pockets formed by the skin and stiffeners shown in figure 17 are isosceles triangles. The 
pockets shown in figure 18 for the Kagome grid are isosceles triangles and hexagons. The 
notation used for the material and section properties and orientation angle of each stiffener family 
is given in Table 3 for both stiffener arrangements. In particular, the subscripts and superscripts 
"s" refer to quantities associated with the stringers, and the subscripts and superscripts "dl" and 
"d2" refer to quantities associated with the two nonidentical families of diagonals. The stiffness 
expressions presented subsequently are obtained from equations (23) and (25)-(28) by applying 
equations (25)-(28) to each family of stiffeners with the attributes given in Table 3. In the 
expressions that follow, the stiffener extensional modulus, shear modulus, eccentricity, moment 
of inertia, torsion constant, and transverse-shear correction factors refer to the corresponding 
effective quantities defined in Appendix A for a nonhomogeneous, specially orthotropic beam. In 

addition, the length of the diagonal stiffeners is given by l = a/l 2 x + L 2 ' • 
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D 2 6 — D 26 + y 3 E dl I dl — E d2 I d2 + I G d2 J d2 — G dl J dl 1 — x 2 

OT ^ /IT 
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where the inplane-shear-deformation parameters are defined as 
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and the transverse-shear-deformation parameters are defined as 
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Setting these shear-deformation parameters equal to zero eliminates the effects of stiffener shear 
deformation. 
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For the special case in which the triangle-shaped pockets enclosed by the stiffeners form 
equilateral triangles and the properties of all stiffeners are identical, L x = L and L y = and 

the stiffnesses reduce to 
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Appendix D 

Beam-Member Stiffness Coefficients 

In the present study, it was found to be convenient to express equation (A25) as 

EsA s 0 0 E s A s y s E s A s z s 0 

• XyE s A s 0 0 0 -XyE s A s z s 

• • x z E s A s 0 0 x z E s A s y s 

• * * Egl^ E s I yz 0 

• * * * E S I YY 0 

symmetric • • • • G S J S 

such that equation (36) yields the following nonzero elements of the symmetric matrix \{? P ] : 


& P U = E S A S cos^slcos 21 ?;; + x Y smV s j (D2) 

<^i 2 = E S A S sin 'Ps cos 2 W s ( 1 - x®) (D3) 

= E S A S sin*, COS2V.) (D4, 

0 P U — E S A S cos 2 l P s ^z s coshp s + XyZ s sin W s j (D5) 

e'l = E S A S sinV s cosV s (z s - x Y f s ) (D6) 

= E S A S sinW s cos 1 ^! cosAb s - ytyZ s cos 2 l P s j (D7) 

^22 = E S A S skf x P s |sin~ l I , s + Xy cos^ j (D8) 

^ = EsA - ^ cos ' I ' ! ( sinV! + ^ cos2>p >) <D9) 
^(4 = E s A s sin 2 ^ coshF s (z s - Xyl s ) (DIO) 
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eh = E S A S sin W s ^z s sin'Eg + r s Y z s cos V s j (D1 1) 

eh = E s A s sin'Fs cosW s |z s sin X F S + ^-TyZ s cos2W s j (D12) 

el, = E S A JsinV, cos 2 W s + ixYCOs 2 2W s ) (D13) 

^ = E > As sin ' I '° cos ' 1 ’* ( Zs cosX “ 5* cos2,p -) <D14) 

els = E S A S sin l P s cos l I f s |z s sin W s + ^-x yf s cos2 v E s j (D15) 

el 6 = E s A s |z s sin"W s cos 2 W s + ^-XyZ s cos 2 2W s j (D16) 

e' 44 = cosV s Je s Iyy cos~W s + G S J S sin ^ s j (D17) 

el 5 = sinV s cosX(e s Iyy - G S J S ) (D18) 

^ l .‘" 1 < di 9 ) 

el 7 = sin 2 W s cosW s E s A s Xz y s (D20) 

els = cos 2 'P s sin'P s E s A s Tz y s (D21) 

el 5 = sinV s (cos 2 W s G S J S + sinV s E s Iy Y ] (D22) 

* = cosW ‘( E ^ sinV > + i > “H <D23) 

eh = - sin~»Eg cos^E s E s A s Xz y s (D24) 

eh = - cos T F s sin^Pg EgAgXz y s (D25) 

eh = E s Iyy SinVg COsVg + ^GgJg COS^Wg (D26) 
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e'h = - -ysinW s cos2 t I , s E s A s x z y s 

(D27) 

e P bS = - icosW s cos2 1 E s E s A s x z y s 

(D28) 

e”-,-, = sin 'Ps E s A s x z 

(D29) 

= sin'Ps cos'Pg E s A s Xz 

(D30) 

^88 = cos 2 W s E s A s x z 

(D31) 

where 


.A k Y G s 
Y _ E 

E s 

(D32) 

is an inplane-shear-deformation parameter and 


.A _ k z G s 

z “ E 
E s 

(D33) 

is a transverse-shear-deformation parameter. Setting these shear-deformation parameters equal 
to zero eliminates the effects of stiffener shear deformation. 
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Appendix E 

Mathematica Program for the Basic-Cell Energy-Equivalence Method 


The computer program 142 used in the present study to determine the contribution of the stiffeners 
to the equivalent-continuum plate stiffnesses, based on the basic-cell energy-equivalence method 
described herein, is presented subsequently. The input used in this program for the orthogonal 
stiffener arrangement with diagonal braces shown in figure 14 is given in Table 4 and the output 
is presented in Appendix B, where a = 1. IGRID = 1 corresponds to the special case in which the 
stiffening grid has no diagonal braces and the stringers are generally different from the ribs (see 
figure 14a). In contrast, IGRID = 2 corresponds to the general case in which the stiffening grid 
has two nonidentical families of diagonal braces in addition to the stringers and ribs. For IGRID 
= 3, all diagonal braces are identical. Similarly, the input used in the program for the orthogonal 
stiffener arrangement with diagonal braces shown in figure 16 is given in Table 5 and the output 
is presented in Appendix B, where a = 1/2. IGRID = 4 corresponds to the general case in which 
the stiffening grid has two nonidentical families of diagonal braces in addition to the stringers and 
ribs. For IGRID = 5, all diagonal braces are identical. 

The input used in this program for the isosceles-triangle stiffener arrangement shown in figure 
17 is given in Table 6 and the output is presented in Appendix C. IGRID = 6 corresponds to the 
general case in which the stiffening grid has two nonidentical families of diagonal braces that are 
different from the stringers (see figure 17a). IGRID = 7 corresponds to the special case of an 
equilateral-triangle grid with identical braces and stringers. Similarly, the input used in this 
program for the Kagome stiffener arrangement, with two nonidentical families of diagonal braces 
that are different from the stringers, shown in figure 18 (IGRID = 8) is given in Table 7 and the 
output was found to be identical to that presented in Appendix C for the isosceles-triangle 
stiffener arrangement. 

The input used in the program for the hexagon- shaped stiffener arrangement shown in figure 
21 is given in Table 8 and the output is presented in Appendix F. IGRID = 9 corresponds to the 
general case in which one pair of opposite sides have a different length and different properties 
than the remaining sides, and the properties of each family of diagonal members are different. 
IGRID =10 corresponds to the special case of an equilateral hexagon-shaped stiffener 
arrangement with identical side lengths and identical diagonal members. 

The input used in the program for the star-shaped stiffener arrangement shown in figure 23 is 
given in Table 9 and the output is presented in Appendix G. IGRID =11 corresponds to the 
general case in which the star shape is made from identical isosceles triangles with coincident 
centroids, and the two families of diagonal members have different properties. IGRID = 12 
corresponds to the special case in which the star shape is made from identical equilateral triangles 
with coincident centroids, and with identical diagonal members. 

IGRID =13 corresponds to the hexagon-shaped core of the stiffened plate shown in figure 25. 
For this case, the core is made of isosceles hexagons formed by beams with rectangular cross 
sections, and is centered on plane z = 0, as shown in figure 25. The basic cell for this case is 
identical to that shown in figure 21b. Likewise, IGRID = 14 is for the special case in which the 
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core is made of isotropic beams with the same length. The output for these cases is presented in 
Appendix H. 


IGRID = 15 corresponds to the star-shaped core of the stiffened plate shown in figure 27. For 
this case, the core is made of isosceles-star shapes formed by beams with rectangular cross 
sections, and is centered on plane z = 0, as shown in figure 27. The basic cell for this case is 
identical to that shown in figure 23b. Likewise, IGRID = 16 is for the special case in which the 
core is made of isotropic beams with the same length. The output for these cases is presented in 
Appendix J. 

In the computer program listing that follows, the user-defined input is given by the single red 
line that indicates a selection for the IGRID variable. The program is then executed by using the 

® 

"evaluate notebook" Mathematica command. For each of the stiffener arrangements defined by 
IGRID, the (x, y) coordinates of the starting point and ending point of members forming a given 
basic cell are arranged in row vectors, in ascending order. For example, the starting and ending 
points, or nodes, of the members in the basic cell shown in figure 14b correspond to the filled red 
circles and are numbered from 1 to 8. The corresponding row vectors for the coordinates of the 
nodes, specified in ascending order from 1 to 8, are 

x = {-Lx/2, Lx/2, 0, 0, -Lx/2, Lx/2, Lx/2, -Lx/2} 

y = { 0, 0, -Ly/2 , Ly/2 , -Ly/2, Ly/2, -Ly/2, Ly/2} 

Next, the starting and ending node number and the geometric and material properties of each 
beam member forming a given basic cell are specified in a table. For example, the table for the 
basic cell shown in figure 14b is given by 


PROPS={{l,2,Es,As, YBs, Y2Bs, ZBs, Z2Bs, IYs, IZs, IYZs,Gs,TauYs,TauZs, Js,Lx} , 

{3 , 4,Er,Ar, YBr, Y2Br, ZBr, Z2Br, IYr, IZr, IYZr,Gr,TauYr,TauZr, Jr,Ly} , 

{5,6, Edl , Adi , YBdl , Y2Bdl , ZBdl , Z2Bdl , IYdl , IZdl , IYZdl , Gdl , TauYdl , TauZdl , Jdl , Ld} , 

{7,8, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , I YZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , Ld} } 

Each row in this table corresponds to one beam member of the basic cell shown in figure 14b. The 
nomenclature follows that given in Table 4 and Appendix D. 

Upon forming the geometry and material properties of the specified basic cell, the strain energy 
of each member in the basic cell is computed by using equation (37) and then summed and divided 
by the area of the basic cell, as indicated by equation (39a). To obtain this equivalent-continuum 
strain-energy density, cos l F s and sin v F s for each beam member in the basic cell are computed by 
using the nodal coordinates the length of each member. This calculation enables formation of the 
transformation matrix given by equation (34a). Then the beam constitutive matrix given by 
equation (Dl) is formed, the strain-equivalence matrix given by equation (30d) is formed, and the 
calculations indicated by equations (32) and (36) are performed. It is important to note that the 
matrix elements of equation (Dl) are the effective, noncentroidal values. All the calculations are 
completed within a "DO LOOP" that corresponds to the summation symbol in equation (39a). 
Upon forming the equivalent-continuum strain-energy density, the contributions of the stiffeners 
to the equivalent-plate stiffnesses are calculated following equations (43). The code that performs 
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these calculations is given as follows. 


/************************************************************************\ 


LIST OF STIFFENER ARRANGEMENTS *) 

*) 

0=1 Orthogonal Grid of ribs and stringers *) 

2 Orthogonal Grid of ribs and stringers with two *) 

nonidentical braces per bay *) 

3 Orthogonal Grid of ribs and stringers with two *) 

identical braces per bay *) 

4 Orthogonal Grid of ribs and stringers with nonidentical *) 

braces and one brace per bay *) 

5 Orthogonal Grid of ribs and stringers with identical *) 

braces and one brace per bay *) 

6 Isosceles-Triangle Grid with nonidentical diagonals *) 

7 Equilateral-Triangle Grid with identical members *) 

8 KAGOME Grid depicted in figure 18 with *) 

nonidentical diagonals *) 

9 Isosceles-Hexagon Grid with nonidentical diagonals *) 

10 Equilateral-Hexagon Grid with identical diagonals *) 

11 Isosceles-star-cell Grid with nonidentical diagonals *) 

12 Equilateral-star-cell Grid with identical diagonals *) 

13 Isosceles-Hexagon Core made of rectangular-cross-section *) 

beams and with identical diagonals *) 

14 Equilateral-Hexagon Core made of rectangular-cross-section *) 

isotropic beams *) 

15 Isosceles-star-cell Core with identical diagonals *) 

16 Equilateral-star-cell Core with identical isotropic members*) 

) 


/ ************************************************************************ 


I GRID =3; (* SELECT GRID-STIFFENER ARRANGEMENT *) 


If [ IGRID == 1, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 
(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers =2, CellArea = Lx Ly, 

x = {-Lx/2, Lx/2, 0, 0), 

y = { 0, 0, -Ly/2 , Ly / 2 } , 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 
(* PROPERTIES FOR THE BASIC CELL *) 


PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , Lx} , 

{3,4, Er , Ar , YBr , Y2Br , ZBr , Z2Br , IYr,IZr, IYZr , Gr , TauYr , TauZr , Jr ,Ly} } }] ; 

If [IGRID == 2, { 


(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 
(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 


NMembers = 4, CellArea = Lx Ly, 


(* 

(* 


x = {-Lx/2, Lx/2, 0, 0, -Lx/2, Lx/2, 

y = { 0, 0, -Ly/2, Ly/2, -Ly/2, Ly/2, 

DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND 
PROPERTIES FOR THE BASIC CELL 


Lx/2, -Lx/2}, 
-Ly/2, Ly/2}, 

ENDING NODES AND *) 

*) 


PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , Lx} , 

{3,4,Er,Ar, YBr, Y2Br, ZBr, Z2Br, I Yr, IZr, IYZr, Gr, TauYr, TauZr, Jr, Ly}, 

{5,6, Edl , Adi , YBdl , Y2Bdl , ZBdl , Z2Bdl , IYdl ,IZdl,IYZdl, Gdl , TauYdl , TauZdl , Jdl , Ld} , 

{7,8, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , IYZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , Ld} } } ] 


If [IGRID==3, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 
(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=4 ,CellArea=Lx Ly, 

x={-Lx/2, Lx/2, 0, 0, -Lx/2, Lx/2, Lx/2, -Lx/2}, 

y={ 0, 0, -Ly/2, Ly/2, -Ly/2, Ly/2, -Ly/2, Ly/2}, 
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(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js,Lx}, 

{3,4, Er,Ar, YBr, Y2Br , ZBr , Z2Br , IYr , IZr , IYZr , Gr, TauYr , TauZr , Jr,Ly}, 

{ 5 , 6 , Ed , Ad , YBd , Y2Bd , ZBd , Z2Bd , I Yd , IZd , IYZd, Gd , TauYd , TauZd , Jd , Ld} , 

{7,8, Ed, Ad , YBd , Y2Bd , ZBd , Z2Bd , I Yd, IZd, I YZd , Gd , TauYd , TauZd , Jd , Ld} > }] ; 


If [IGRID==4, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=6 ,CellArea=4Lx Ly, 

x={-Lx, Lx, 0, 0, -Lx, Lx, Lx, -Lx}, 
y={ 0, 0, -Ly, Ly, -Ly, Ly, -Ly, Ly} , 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , 2Lx} , 

{ 5 , 7 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , 2Lx} , 

{3,4, Er , Ar , YBr, Y2Br , ZBr , Z2Br , IYr, IZr, IYZr , Gr , TauYr , TauZr , Jr , 2Ly } , 

{5,8, Er , Ar , YBr, Y2Br , ZBr, Z2Br , IYr, IZr, IYZr , Gr , TauYr , TauZr , Jr , 2Ly} , 

{5,6, Edl , Adi , YBdl , Y2Bdl , ZBdl , Z2Bdl , IYdl , IZdl , IYZdl , Gdl , TauYdl , TauZdl , Jdl , 2Ld} , 

{7,8, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , IYZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , 2Ld} } } ] ; 


If [IGRID==5, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=6 ,CellArea=4Lx Ly, 

x={-Lx, Lx, 0, 0, -Lx, Lx, Lx, -Lx}, 
y={ 0, 0, -Ly, Ly, -Ly, Ly, -Ly, Ly}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , 2Lx} , 

{5,7 ,Es, As, YBs, Y2Bs, ZBs, Z2Bs, IYs, IZs, IYZs, Gs, TauYs, TauZs, Js,2Lx} , 
{3,4,Er,Ar, YBr, Y2Br, ZBr, Z2Br, IYr, IZr, IYZr, Gr, TauYr, TauZr, Jr, 2Ly} , 
{5,8,Er,Ar, YBr, Y2Br, ZBr, Z2Br, IYr, IZr, IYZr, Gr, TauYr, TauZr, Jr, 2Ly} , 

{5, 6, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd,2Ld} , 

{ 7 , 8 , Ed , Ad , YBd , Y2Bd , ZBd , Z2Bd , IYd , IZd , I YZd , Gd , TauYd , TauZd , Jd , 2Ld} } } ] ; 


If [IGRID==6, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=3 , CellArea=Lx Ly , 

x={-Lx/2, Lx/2, -Lx/4, Lx/4, Lx/4, -Lx/4}, 

y={ 0, 0, -Ly/2, Ly/2, -Ly/2, Ly/2}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , Lx} , 

{3, 4, Edl, Adi, YBdl, Y2Bdl, ZBdl, Z2Bdl, IYdl, IZdl, IYZdl, Gdl, TauYdl, TauZdl, Jdl, Ld} , 

{5,6, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , I YZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , Ld} } } ] ; 


If [ IGRID==7 , { 

(* SET PROPERTY EQUALITIES FOR IDENTICAL MEMBERS AND SPECIAL CASE CRITERIA *) 

Lx=L , Ly=Sqrt [ 3 ] /2 L , Ld=L , Ad=As , Ed=Es , Gd=Gs , Jd= Js , 

YBd=YBs , Y2Bd=Y2Bs , ZBd=ZBs , Z2Bd=Z2Bs , IYd=IYs , IZd=IZs , IYZd=IYZs , 

TauYd=TauYs , TauZd=TauZs , 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=3 ,CellArea=Lx Ly, 

x={-Lx/2, Lx/2, -Lx/4, Lx/4, Lx/4, -Lx/4}, 

y={ 0, 0, -Ly/2, Ly/2, -Ly/2, Ly/2}, 
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(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js,Lx}, 

{ 3 , 4 , Ed , Ad , YBd , Y2Bd , ZBd , Z2Bd , I Yd , IZd , IYZd, Gd , TauYd , TauZd , Jd , Ld} , 

{5,6, Ed , Ad, YBd, Y2Bd , ZBd , Z2Bd , IYd , IZd , I YZd , Gd , TauYd , TauZd , Jd , Ld} > }] ; 


If [IGRID==8, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=4 ,CellArea=2 Lx Ly, 

x={-Lx/2, Lx/2, -Lx/2, Lx/2, -Lx/2, Lx/2, Lx/2, -Lx/2}, 
y={-Ly/2, -Ly/2, Ly/2, Ly/2, -Ly, Ly, -Ly, Ly) , 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1 , 2 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , Lx} , 

{3 , 4 , Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , Lx} , 

{5,6, Edl , Adi , YBdl , Y2Bdl , ZBdl , Z2Bdl , IYdl ,IZdl,IYZdl, Gdl , TauYd 1 , TauZdl , Jdl , 2L} , 

{7,8, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , IYZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , 2L} } } ] ; 


If [IGRID==9, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=5 ,CellArea=2*a* (b+c) , 

x={ a/2, 0, -a/2, a/2, 0, -a/2}, 

y={ - (b+c) /2 , -c/2, - (b+c) /2 , (b+c)/2, c/2, (b+c)/2}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ {2 , 1 , Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , IYZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , L/2 } , 

{2,3, Edl , Adi , YBdl , Y2Bdl , ZBdl , Z2Bdl , IYdl ,IZdl,IYZdl, Gdl , TauYd 1 , TauZdl , Jdl , L/2 } , 

{2,5, Er , Ar , YBr , Y2Br , ZBr , Z2Br , IYr,IZr, IYZr , Gr , TauYr , TauZr , Jr , c} , 

{5, 4, Edl, Adi, YBdl, Y2Bdl, ZBdl, Z2Bdl, IYdl, IZdl, I YZdl, Gdl, TauYdl, TauZdl, Jdl, L/2}, 

{5,6, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , I YZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , L/2 } } } ] 

If [IGRID==10, { 

(* SET PROPERTY EQUALITIES FOR IDENTICAL MEMBERS AND SPECIAL CASE CRITERIA *) 
a=Sqrt[3]/2 L,b=L/2,c=L, 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=5 ,CellArea=2*a* (b+c) , 

x={ a/2, 0, -a/2, a/2, 0, -a/2}, 

y={ - (b+c) /2 , -c/2, - (b+c) /2 , (b+c)/2, c/2, (b+c)/2}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ {2, 1, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, L/2} , 

{ 2 , 3 , Ed , Ad , YBd , Y2Bd , ZBd , Z2Bd , IYd , IZd , I YZd , Gd , TauYd , TauZd , Jd , L/2 } , 

{2,5,Er,Ar, YBr, Y2Br, ZBr, Z2Br, IYr, IZr, IYZr, Gr, TauYr, TauZr, Jr, c}, 

{5, 4, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, L/2} , 

{ 5 , 6 , Ed , Ad , YBd , Y2Bd , ZBd , Z2Bd , IYd , IZd , I YZd , Gd , TauYd , TauZd , Jd , L/2 } } } ] ; 

If [IGRID==11, { 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=12 , CellArea=4 B H/3, 

x={B/3, B/2, B/6, 0, -B/6, -B/2, -B/3, -B/2, -B/6, 0, B/6, B/2}, 

y={ 0, H/3, H/3, 2 H/3, H/3, H/3, 0, -H/3, -H/3, -2 H/3, -H/3, -H/3}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1,12, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , IYd2 , IZd2 , IYZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , Ld/3 } , 
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{1,2, Edl , Adi , YBdl , Y2Bdl ,ZBdl,Z2Bdl, IYdl, IZdl, IYZdl, Gdl , TauYdl , TauZdl , Jdl , Ld/3 } , 
{3,2, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , B/3 } , 

{3,4, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , 1 Yd2 , IZd2 , IYZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , Ld/ 3 } , 
{5, 4, Edl, Adi, YBdl, Y2Bdl , ZBdl , Z2Bdl , IYdl , IZdl , IYZdl , Gdl, TauYdl, TauZdl, Jdl, Ld/3}, 
{5,6, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3) , 

{7,6, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , 1 Yd2 , IZd2 , IYZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , Ld/ 3 } , 
{7, 8, Edl, Adi, YBdl, Y2Bdl, ZBdl, Z2Bdl, IYdl, IZdl, IYZdl, Gdl, TauYdl, TauZdl, Jdl, Ld/3}, 
{9,8, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} , 

{9,10, Ed2 , Ad2 , YBd2 , Y2Bd2 , ZBd2 , Z2Bd2 , 1 Yd2 , IZd2 , 1 YZd2 , Gd2 , TauYd2 , TauZd2 , Jd2 , Ld/ 3 } 
{11, 10, Edl , Adi , YBdl , Y2Bdl , ZBdl , Z2Bdl , IYdl, IZdl, IYZdl, Gdl , TauYdl , TauZdl , Jdl , Ld/3 } 
{ll,12,Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} } }] ; 

If [IGRID==12,{ 

H=Sqrt [ 3 ] L/2,B= L,Ld=L, 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND * ) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=12 , CellArea=4 B H/3, 

x={B/3 , B/2 , B/6, 0, -B/6, -B/2, -B/3, -B/2, -B/6, 0, B/6, B/2}, 

y={ 0, H/3, H/3, 2 H/3, H/3, H/3, 0, -H/3, -H/3, -2 H/3, -H/3, -H/3}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

PROPS={ { 1,12, Ed, Ad, YBd, Y2Bd, ZBd , Z2Bd , I Yd , IZd, IYZd, Gd, TauYd , TauZd, Jd, Ld/3} , 

{1,2, Ed, Ad , YBd , Y2Bd , ZBd , Z2Bd , IYd , IZd , IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{3,2, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} , 

{3 ,4, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, I YZd,Gd, TauYd, TauZd, Jd,Ld/3} , 

{5,4, Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{5,6, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} , 

{7 ,6, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, I YZd,Gd, TauYd, TauZd, Jd,Ld/3} , 

{7,8, Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{9,8, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} , 

{9,10, Ed, Ad, YBd, Y2Bd,ZBd,Z2Bd,IYd,IZd, IYZd, Gd, TauYd , TauZd, Jd, Ld/3}, 

{11, 10, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3} , 

{ll,12,Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} } }] ; 

If [IGRID==13,{ 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS * ) 

NMembers=5 ,CellArea=2*a* (b+c) , 

x={ a/2, 0, -a/2, a/2, 0, -a/2}, 

y={ - (b+c) /2 , -c/2, -(b+c)/2, (b+c)/2, c/2, (b+c)/2}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

Ad=hc td, 

Ar=hc tr, 

ZBd=0 , 

Z2Bd=0 , 

ZBr=0 , 

Z2Br=0 , 

IYd=hc'3 td/12 , 

IZd=td'3 hc/12 , 

IYZd=0 , 

IYr=hc'3 tr/12 , 

IZr=tr A 3 hc/12, 

IYZr=0 , 

PROPS={{2, 1, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd,L/2} , 

{2,3 , Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd,L/2}, 

{2,5, Er , Ar , YBr , Y2Br , ZBr , Z2Br , IYr , IZr , IYZr , Gr, TauYr, TauZr, Jr , c} , 

{5, 4, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd,L/2}, 

{5, 6, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd,L/2} } }] ; 

If [IGRID==14,{ 

a=Sqrt [ 3 ] / 2 L,b=L/2,c=L, 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 


NMembers=5 ,CellArea=2*a* (b+c) , 
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x={ a/2, 0, -a/2, a/2, 0, -a/2}, 

y={ - (b+c) /2 , -c/2, -(b+c)/2, (b+c)/2, c/2, (b+c)/2}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

Ad=hc td, 

Ar=hc tr, 

ZBd=0 , 

Z2Bd=0 , 

ZBr=0 , 

Z2Br=0 , 

IYd=hc~3 td/ 12 , 

IZd=td'3 hc/12 , 

IYZd=0 , 

IYr=hc~3 tr/12 , 

IZr=tr'3 hc/12, 

IYZr=0 , 

Ed=EE , Er=EE , Gd=G, Gr=G, TauYd=k G/EE, TauYr=k G/EE , TauZd=k G/EE , TauZr=k G/EE, 
Jd=td~3 hc/3 (1-Kd) , Jr=tr~3 hc/3 (1-Kr), 

PROPS={ {2 , 1 , Ed, Ad , YBd , Y2Bd , ZBd , Z2Bd , IYd , IZd , IYZd, Gd, TauYd, TauZd, Jd, L/2 } , 

{2,3, Ed, Ad , YBd , Y2Bd , ZBd , Z2Bd , lYd , IZd , IYZd, Gd, TauYd, TauZd, Jd, L/2 } , 

{2,5, Er , Ar , YBr , Y2Br , ZBr , Z2Br , IYr , IZr , IYZr , Gr, TauYr, TauZr, Jr , c} , 

{5,4, Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, L/2 } , 

{5,6, Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, L/2 } } }] ; 

If [IGRID==15,{ 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND * ) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS * ) 

NMembers=12 , CellArea=4 B H/3, 

x={B/3, B/2 , B/6, 0, -B/6, -B/2, -B/3, -B/2, -B/6, 0, B/6, B/2}, 

y={ 0, H/3, H/3, 2 H/3, H/3, H/3, 0, -H/3, -H/3, -2 H/3, -H/3, -H/3}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 

As = ts he, Ad = td hc,IYs = ts hc A 3/12,IYd = td hc~3/12, 

ZBd=0 , ZBs=0 , Z2Bd=0 , Z2Bs=0, Js = ts~3 hc/3 (l-Ks),Jd = td"3 hc/3 (1-Kd), 

PROPS={ { 1,12, Ed , Ad, YBd, Y2Bd , ZBd , Z2Bd , IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3} , 

{1,2, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd,Ld/3} , 

{3,2, Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , B/3} , 

{3,4, Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{5, 4, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3} , 

{5,6, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} , 

{7,6, Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{7 ,8, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, I YZd,Gd, TauYd, TauZd, Jd,Ld/3} , 

{9,8, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} , 

{9,10, Ed , Ad, YBd, Y2Bd, ZBd, Z2Bd , IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3}, 

{11, 10, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3} , 
{ll,12,Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} } }] 

If [IGRID==16,{ 


B=L , H=Sqrt [ 3 ] L/2, 

(* DEFINE THE NUMBER OF MEMBERS, THE AREA OF THE BASIC CELL, AND *) 

(* THE SEQUENTIAL GLOBAL (x, y) NODAL COORDINATES AS ROW VECTORS *) 

NMembers=12 , CellArea=4 B H/3, 

x={B/3, B/2, B/6, 0, -B/6, -B/2, -B/3, -B/2, -B/6, 0, B/6, B/2}, 

y={ 0, H/3, H/3, 2 H/3, H/3, H/3, 0, -H/3, -H/3, -2 H/3, -H/3 ,-H/3}, 

(* DEFINE A TABLE OF BEAM-MEMBER BEGINNING AND ENDING NODES AND *) 

(* PROPERTIES FOR THE BASIC CELL *) 


Ld=L, 

As = t he, Ad = t he, IYs = t hc~3/12,IYd = t hc"3/12, 

ZBd=0 , ZBs=0 , Z2Bd=0 , Z2Bs=0, Js = t~3 hc/3 (l-KK),Jd = t"3 hc/3 (1-KK), 

Es=EE , Ed=EE , Gs=G, Gd=G, TauYd=k G/EE, TauYs=k G/EE , TauZd=k G/EE,TauZs=k G/EE, 

PROPS={ { 1,12, Ed , Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd , IYZd, Gd, TauYd, TauZd, Jd, Ld/3} , 
{1,2, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3} , 
{3,2, Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} , 
{ 3 , 4 , Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, IYd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 
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{5, 4, Ed, Ad, YBd, Y2Bd, ZBd, Z2Bd, I Yd, IZd, I YZd,Gd, TauYd, TauZd, Jd,Ld/3> , 

{5,6, Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs , IZs , IYZs , Gs , TauYs , TauZs , Js , B/3 } , 

{7,6, Ed, Ad , YBd , Y2Bd , ZBd, Z2Bd, I Yd , IZd , IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{7,8, Ed, Ad , YBd, Y2Bd, ZBd, Z2Bd, I Yd , IZd , IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{9,8, Es , As , YBs , Y2Bs , ZBs , Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3) , 

{9,10, Ed , Ad, YBd, Y2Bd , ZBd, Z2Bd, I Yd, IZd, IYZd, Gd, TauYd, TauZd, Jd, Ld/3}, 
{11,10, Ed , Ad, YBd, Y2Bd , ZBd , Z2Bd , I Yd , IZd , IYZd, Gd, TauYd, TauZd, Jd, Ld/3 } , 

{11 , 12 , Es , As , YBs ,Y2Bs,ZBs, Z2Bs , IYs, IZs, IYZs, Gs , TauYs , TauZs , Js , B/3} } }] 


^************************************************************************J 

(* BEGIN A DO-LOOP THAT COMPUTES THE STRAIN ENERGY CONTRIBUTION OF EACH *) 
(* BEAM MEMBER TO THE STRAIN ENERGY OF THE EQUIVALENT CONTINUUM *) 

( * *) 
(* THE CURRENT BEAM MEMBER IN THE LOOP IS IDENTIED BY THE INDEX K *) 

^************************************************************************j 

(* INITIALIZE THE STRAIN-ENERGY OF THE EQUIVALENT CONTINUUM *) 

StrainEnergy =0; 

Do [ 


(* COMPUTE GLOBAL COORDINATES OF THE CURRENT BEAM MEMBER IJ *) 


xl = x[[ 
XJ = x[[ 
Yl = Y[[ 
y J = y [ [ 


PROPS [[K, 1]] ]]; 
PROPS [ [ K , 2]] ]]; 
PROPS [[K, 1]] ]]; 
PROPS [ [ K , 2]] ]]; 


(* COMPUTE THE LENGTH OF THE CURRENT BEAM MEMBER *) 


LengthlJ = PROPS[[K, 16]]; 

(* COMPUTE m = COS(psi) and n = SIN(psi) FOR THE CURRENT BEAM MEMBER *) 

m = (xJ - xl ) /LengthlJ; 
n = (yJ - yl ) /LengthlJ; 


(* COMPUTE THE PROPERTIES OF THE CURRENT BEAM MEMBER *) 


EX=PROPS [ [K,3] ] ; 

A=PROPS [ [K,4] ] ; 

Ybar=PROPS [ [K, 5] ] ; 

Y2bar=PROPS[ [K,6] ] ; 

Zbar=PROPS [ [K,7] ] ; 

Z2bar=PROPS[ [K,8] ] ; 

IYY=PROPS [ [K,9] ] ; 

IZZ=PROPS [ [K, 10] ] ; 

IYZ=PROPS [ [K, 11] ] ; 

G=PROPS [ [K,12] ] ; 

TauY=PROPS [ [K, 13] ] ; 

TauZ=PROPS [ [K, 14] ] ; 

J=PROPS [ [K,15] ] ; 

(* COMPUTE THE STIFFNESS MATRIX FOR THE CURRENT BEAM MEMBER; see equation (Dl) *) 


Cll 

= 

EX 

A; 

C14 

= 

EX 

A Ybar; 

C15 

= 

EX 

A Zbar; 

C22 

= 

EX 

A TauY; 

C26 

= 

-EX A TauY Z2bar 

C33 

= 

EX 

A TauZ; 

C36 

= 

EX 

A TauZ Y2bar; 

C44 

= 

EX 

IZZ; 

C45 

= 

EX 

IYZ ; 

C55 

= 

EX 

IYY ; 

C66 

= 

G v 

J; 



Cll, 

0 , 

o. 

C14 , 

C15 , 

0 }, 


0 , 

C22 , 

o. 

o. 

o. 

C26 } , 


0 , 

o. 

C33 , 

o. 

0 , 

C36} , 


C14 , 

0 , 

o. 

C44 , 

C45 , 

0 }, 

{ 

C15 , 

o. 

o. 

C45 , 

C55 , 

o>. 


0 , 

C26 , 

C36, 

o. 

o. 

C66} } 
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(* COMPUTE THE TRANSFORMED STIFFNESS MATRIX FOR THE CURRENT BEAM MEMBER; see equations (30) -(36) *) 


Ematrix = 

{ { 

1 , 0 , 

o. 

o 

o 

O 

<3 

<3 





{ 

o 

o 

1/2, 

o 

o 

0 , 0 , 0>, 





{ 

o 

o 

0 / 

o 

o 

0 , 0 , 1>, 





{ 

o 

o 

o. 

o 

o 

0 , 0 , 0>, 





{ 

o 

o 

o. 

1 , 0 , 

0 , 0 , 0>, 





{ 

o 

o 

o. 

i 

o 

o 

1/2, 0, 0} >; 




CPmatrix = 

Transpose [Ematrix] .Cmatrix 

.Ematrix; 





Transformationmatrix = 

{ 







{ m A 2 , 

n~2. 


m*n. 

0 f 

0 , 

0 , 

0 , 

0}, 

{ n~2, 

m~2 , 


m*n. 

o. 

o. 

0 , 

0 , 

0>, 

{ -2*m*n, 

2*m*n, 

m~ 

2 - n~2. 

o. 

0 , 

0 , 

0 , 

0>, 

{ 

0 , 

0 , 


o. 

m^2 , 

n~2. 

m*n. 

0 , 

0}, 

{ 

0 , 

0 , 


0 

n^2 , 

m~ 2 , 

m*n. 

0 , 

0>, 

{ 

0 , 

0 , 


o. 

-2*m*n, 

2*m*n, 

m~2 - n A 2, 

0 , 

0}, 

{ 

0 , 

0 , 


o. 

o. 

0 , 

0 , 

m. 

-n > , 

{ 

0 , 

0 , 


0, 

0 , 

0 , 

0 , 

n. 

m> > 


Cpmatrix = Transpose [Transf ormationmatrix] . CPmatrix . Transformationmatrix; 


(* COMPUTE STRAIN ENERGY CONTRIBUTION OF THE CURRENT BEAM MEMBER; see equation (37) *) 
StrainList = {epxx, epyy, gamxy, kapxx, kapyy, kapxy, gamyz, gamxz}; 

MemberStrainEnergy = LengthIJ/2 StrainList .Cpmatrix. StrainList; 

(* ADD THE STRAIN-ENERGY CONTRIBUTION OF THE CURRENT BEAM MEMBER TO THE *) 

(* CONTINUUM STRAIN ENERGY; see equation (39a) *) 

StrainEnergy = StrainEnergy + MemberStrainEnergy, 

(* CLOSE DO-LOOP *) {K, 1, NMembers}] 

(* COMPUTE THE TOTAL STRAIN ENERGY DENSITY OF THE EQUIVALENT CONTINNUM; see equation (39a) *) 
StrainEnergyDensity = StrainEnergy/CellArea; 


(* COMPUTE THE A-MATRIX OF FIRST-ORDER TRANSVERSE-SHEAR DEFORMATION PLATE THEORY; see equation (43a) *) 


Print [ 


Astiffenerll = ", FullSimplify [Factor [ D [ StrainEnergyDensity , epxx, epxx] ] ] ]; 


Print [ 

' ”] 

Print [ " 




Print [ 

'Astif fenerl2 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, epxx. 

epyy] ] ] ] ; 

Print [ 

' ”] 

Print [ " 

"]? 



Print [ 

'Astif fener 16 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, epxx. 

gamxy] ] ] ]; 

Print [ 

' ”] 

Print [ " 




Print [ 

'Astif fener 2 2 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, epyy. 

epyy] ] ] ] ; 

Print [ 

' "i 

Print [ " 

"]? 



Print [ 

'Astif fener26 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, epyy. 

gamxy] ] ] ]; 

Print [ 

' "i 

Print [ " 

"]; 



Print [ 

'Astif fener 6 6 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, gamxy 

gamxy ] ] ] ] ; 

Print [ 

' "i 

Print [ " 

"]; Print [ " "]; 



(* COMPUTE 

THE B-MATRIX OF FIRST-ORDER TRANSVERSE-SHEAR DEFORMATION PLATE THEORY; see equation 

Print [ 

'Bstif fenerll 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, epxx. 

kapxx] ] ] ]; 

Print [ 

' "i 

Print [ " 

"]? 



Print [ 

'Bstif fenerl2 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, epxx. 

kapyy] 1 1 ]; 

Print [ 

' "i 

Print [ " 

"]? 



Print [ 

'Bstif fenerl6 

= ", FullSimplify [Factor [ 

D[ StrainEnergyDensity, epxx. 

kapxy ] ] ] ]; 

Print [ 

' "i 

Print [ " 

"]? 
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Print [ 
Print [ 

'Bstif fener22 
' "]; Print [" 


, FullSimplify [Factor [ 

D [ StrainEnergyDensity , 

e PYY / 

kapyy] 

Print [ 
Print [ 

'Bstif fener26 
' "]; Print [" 

"] ; 

, FullSimplify [Factor [ 

D [ StrainEnergyDensity , 

ePYY/ 

kapxy] 

Print [ 
Print [ 

'Bstif fener66 
' " ]; Print [" 


, FullSimplify [Factor [ 
Print [ " "]; 

D [ StrainEnergyDensity , 

gamxy 

kapxy] 


(* COMPUTE THE D-MATRIX OF FIRST-ORDER TRANSVERSE-SHEAR DEFORMATION PLATE THEORY; see equation (43c) *) 


Print [ 
Print [ 

'Dstif fenerll = " 
' " ]; Print [" " ] ; 

, FullSimplify [Factor [ 

D [ StrainEnergyDensity , 

kapxx. 

kapxx] 

Print [ 
Print [ 

'Dstif fenerl2 = " 
' " ]; Print [" " ] ; 

, FullSimplify [Factor [ 

D [ StrainEnergyDensity , 

kapxx. 

kapyy ] 

Print [ 
Print [ 

'Dstif fenerl6 = " 
' "]; Print [" " ] ; 

, FullSimplify [Factor [ 

D [ StrainEnergyDensity , 

kapxx. 

kapxy] 

Print [ 
Print [ 

'Dstif fener22 = " 
' " ]; Print [" " ] ; 

, FullSimplify [Factor [ 

D [ StrainEnergyDensity , 

kapyy. 

kapyy] 

Print [ 
Print [ 

'Dstif fener26 = " 
' " ]; Print [" " ] ; 

, FullSimplify [Factor [ 

D [ StrainEnergyDensity , 

kapyy. 

kapxy] 

Print [ 
Print [ 

'Dstif fener66 = " 
' " ]; Print [" " ] ; 

, FullSimplify [Factor [ 
Print [ " "]; 

D [ StrainEnergyDensity , 

kapxy. 

kapxy] 


(* COMPUTE THE TRANSVERSE -SHEAR STIFFNESS MATRIX OF FIRST-ORDER TRANSVERSE -SHEAR DEFORMATION PLATE 
THEORY; see equation (43d) *) 


Print [ 
Print [ 

'Astif fener44 = " 
' " ]; Print [" " ] ; 

FullSimplify [Factor [ 

D [ StrainEnergyDensity, 

gamyz , 

gamyz ] 

Print [ 
Print [ 

'Astif fener45 = " 
' "]; Print [" " ] ; 

FullSimplify [Factor [ 

D [ StrainEnergyDensity, 

gamyz , 

gamxz] 

Print [ 

'Astif fener55 = " 

FullSimplify [Factor [ 

D [ StrainEnergyDensity, 

gamxz , 

gamxz] 


Print [ " "]; Print[" "]; Print [" " ] ; 

(* RE- INITIALIZE MATHEMATICA *) 

Remove [ " Global ' * " ] ; 

(* CHANGE INPUT AND SELECT "EVALUATE NOTEBOOK" UNDER THE EVALUATION MENU ON THE MATHEMATICA WINDOW HEADER 
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Appendix F 

Equivalent-Plate Stiffnesses for a Plate Reinforced with a Hexagon-Shaped 

Stiffener Arrangement 

The expressions presented in this appendix are for a general laminated-composite plate that is 
stiffened with ribs and two nonidentical families of diagonal braces, as shown in figure 21. The 
ribs and diagonals are eccentric with respect to the plate midplane and the pockets formed by the 
skin and stiffeners shown in figure 21 are isosceles hexagons; that is, the two pairs of diagonal 
members have the same length, which is generally different from the length of the two identical 
ribs. The notation used for the material and section properties and the orientation angle of each 
stiffener family is given in Table 8. In this table and the expressions that follow, the subscript and 
superscript "r" refers to the ribs and the subscripts and superscripts "dl " and "d2" refer to the two 
nonidentical families of diagonals. The stiffness expressions presented subsequently are obtained 
from equations (23) and (25)-(28) by applying equations (25)-(28) to each family of stiffeners 
with the attributes given in Table 8. In the expressions that follow, the stiffener extensional 
modulus, shear modulus, eccentricity, moment of inertia, torsion constant, and the transverse- 
shear correction factors refer to the corresponding effective quantities defined in Appendix A for 
a nonhomogeneous, specially orthotropic beam. 


plate 3 

A n =A n + — 

2(b + c)L 



( 2 A 


f 2 A 


E d iA dl j 

i+4‘A 

| + E d2 A d2 

1+t' v 3 ^ 



l a ' J 


l a ' ) 



a -A P ‘ a,e + ab 

12 — /A 12 ' 


2 b + c L L 


|A dl ( 1 - Ty ) + E d2 A d2 ( 1 - Ty) 


plate 3 b 

A, 6 = A I6 +- —3 

2(b + clL 



!■ 

dl 

2 

\ 

( 62 

2 

-r 


E d iA dl 

+ T 

b — l 
2 1 
a 

1 — E d2 A d2 | 


b 

2 

a 

- 1 



V 






k 


A - A P ‘ ate i EfArC i b 

i v 22 i v 22 ■ z \ i 


4 4 P h,e b 

A 26 — A 26 + — 3 


2a(b + c) 2a(b + c)L 

/ 


, 2 1 


E dl Aj 1 + t 


2(b + c)L~ 


E dl A dl 


1 + f 


Y b 2 


+ E rt9 A., 


< a 2 \ 

9 d2 a 

1 + X Y 

b 


- 1 


-FA 

L -‘A2 n - d2 


1 X i2 

1 + T 


- 1 


(FI) 

(F2) 

(F3) 

(F4) 

(F5) 
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plate E r A r X Y c 

A« + 7 r + 


8a(b 


+ c 


ab” 

2(b + c)L 3 


E d iA dl 


l+ T 


r d 7 2 u 2 A 

x y a - b 


ab 


+ e H9 a (1 


i+ tI 


' 2 ,2' 

a - b 

ab 


(F6) 


D _ T3 plate , 3 

E> ll — -t> 1 1 + _ 3 


2(b + c)L 



7 

2 

A 


f 1 2 A 


E d iA dl 

f _ 

z dl 

dl= b 
+ Xy Z dl — 2 


| + E d2 A d2 j 

d2= b 

Z d2 + X Y Z d2 — 2 



i 

a 

i 


\ 3 / 



(F7) 


D _ i 3 pla,e , ab 

E> 12 — 12 + __ 3 


2(b + c)L L 


lA dl Z d ,| 1 Ty Z dl j + E d2 A d2 Z d2 | 1 'CyZjaj 


(F8) 


_ _ plate a b 

Bi6 = B 16 + 


2(b + c)L 



{ 

dl= 

\ 2 




/ 

d2= 

\ 2 


A 

E d] A dl 

Z dl 

| ^Y Z dl 
2 

b 

2 

a 

- 1 


EdaA d2 

Z d2 

| ^Y Z d2 
2 

b 

2 

a 

- 1 



V 




/ 


V 




7 


(F9) 


plate E,A r CZ r 
D 22 — D 22 "■ 7 r + 


2a(b 


+ c 


2a(b + c)L~ 



7 

2 

\ 


7 2 \ 


E d iA dl j 

IN 

d i= a 

E X Y Z dl 2 


| + E d2 A d2 j 

d2= a 

Z d2 E Xy Z d2 — j 



l 

b 

J 


i b J 



(F10) 


13 _ D P ‘ a,e , b 

b>26 — b> 26 + _ 3 


2(b + c)L 



7 

dl= 

2 


A 


7 

d2= 

2 


A 

E d iA dl 

IN 

| X Y Z dl 

a 

2 

- 1 


E^A d2 

IN 

| X Y Z d2 

a 

2 

- 1 



l 

2 

b 


J 


l 

2 

b 


7 


(FI 1) 


plate E,A r C r = 
B d 6 _ B 66 + , XyZ r 4“ 

8a(b + c) 


ab” 


2(b + c)L” 


3 \ E dl A dl 


z,„ + 


dl= 7 2 i 2 ' 

x y Z tii a — b 
4 ab 


4" E d2 A d2 


Z d ? + 


d2= f 2 .2 

T-y Z d2 a — D 

4 ab 


(FI 2) 


_ _ plate a 

Du = D„ + — 3 


2(b + c)L” 


EdiFi E E d2 I d2 + I GdiJdi + G d2 J d2 J 2 
v ' a 


(FI 3) 
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D 12 -D 12 + 3 E dl I dl + E d 2 I d2 G dl J dl G d 2 J d2 

2(b + c)L L 


D, 6 -D i6 + 


plate 3 b 


2(b + c)L 


3 E dl I dl E d 2 I d2 + 2 


1 fh 2 ) 

9 (G dl J dl — G d 2 J d2 ] ~2 - 1 

2V 'l a i 


_ plate E I C K /_ \ 3 

D 22 = D 22 + . y ^ 7 7 3 E dl I dl + E d2 I d2 + G dl J dl + G d2 J d2 l— 2 

2a(b + c) 2a(b + c)L"[ V 'b 


D -D PlaK + b 
J - y 26 — -^26 ' ( 


2(b + c)L 


E d il d i E d 2 I d2 + 2 (G d ,J dl G d2 J 


b 2 


plate c G,J r t ab 

U 66 - U 66 + ~ 77 . \ + 


8a(b + c) 2(b + c)L 


3 E dl I dl + E d2 I d2 + |G dl J dl + G d2 J d2 


2 u 2 

a - b 


A — A plate 

-/v 44 — 44 H“ 


(E dl A dl Xz + E d2 A d2 x z 2 )b' + EA r x; cL 
2a(b + c)L 


» » p* ate 

A 45 — A 45 + 


E dl A dl x z - E d 2 A d 2 x z 2 b 


2(b + c)L 


» » p ,ate 

A 55 — A 55 + 


-m — , A dl i * d2 

E d iA dl x z + E d 2 A d 2 x z a 


2(b + c)L 


where the inplane-shear-deformation parameters are defined as 


r k Y G r 


di _ k Y G dl 


(F23) 


r d2 _ k Y G d2 


(F24) 
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and the transverse-shear-deformation parameters are defined as 


r k z G r 


Tjy = — 

E r 

(F25) 

dl 

^di _ k z G dl 

z “ E 

-* — ' d 1 

(F26) 

d2 

^d2 _ k z G d2 

z - E 

*^d2 

(F27) 


Setting these shear-deformation parameters equal to zero eliminates the effects of stiffener shear 
deformation. 


For the special case in which all stiffeners have the same length, a = A L , b = L/2, and c = L. 
For this simplification and for identical families of diagonals, the stiffnesses reduce to 


A„ = A n + 


plate vAE d A d [ 1 t y 




4L 


(F22) 


A „ = A'“ + 1 - t{) 


(F23) 


A A P late 

A l6 = A 16 


A — A plate 
za. 22 — za 22 I 


V3" 

36L 


8E r A r + E d A d (l + 3 t y ) 


(F24) 

(F25) 


A - A P ‘ a,e 
^ 26 — ^ 26 


(F26) 


plate t/J 

a 66 = a 66 + 


12L 


2E r A r Ty ^ A , Ty 

+ E d A d 1 H 

3 3 


(F27) 


B n - B„ + 


plate V3E d A d / = _ Xyf d 


4L | + 


(F28) 


D _ -r> pla,e . ^EjAj/ d = \ 

®i2 — B 12 + 12L ' d TyZd J 


(F29) 
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b 16 = b; 


(F30) 


_ nl"" , 7L 

22 22 36L 


8E r ArZ r + E d A d z d (z d + 3x7) 


R _ o plate 

^26 — ^26 


plate VJ 2 EA r XyZ r , B A , TyZ d 

B “- B “ + T2L 3 + E ' AJz ‘ + ^ 


D 11 = Dr + mi j + ^ 


^ plate V3 / \ 

Dl2 = Dl2 +— (E d l d -G d j d ) 


D 16 = d7 


plate V3" 

D22 = e>22 + m: 


'sef + e^ + sga) 


plate 

D26 = D26 


D« = or + ^-(3E d I d + 2G Jr + G d J d ) 


A 44 = A 7 + ^(E d A d Xz + 2E r A r Xz 


A - A P ‘ a,e 

-^-45 -^-45 


plate V^3 E d A d X z 

A„ = A„ 


where the subscript and superscript "d" indicates the properties of the diagonal stiffeners. 


129 



Appendix G 

Equivalent-Plate Stiffnesses for a Plate Reinforced with a Star-Cell-Shaped 

Stiffener Arrangement 

The expressions presented in this appendix are for a general laminated-composite plate that is 
stiffened with stringers and two nonidentical families of diagonal braces, as shown in figure 23. 
The stringers and diagonals are eccentric with respect to the plate midplane and the pockets 
formed by the skin and stiffeners shown in figure 23 are star shapes, formed by isosceles triangles, 
and isosceles hexagons. The two pairs of diagonal members have the same length, which is 
generally different from the length of the two identical stringers. The notation used for the 
material and section properties and orientation angle of each stiffener family is given in Table 9. 
The stiffness expressions presented subsequently are obtained from equations (23) and (25)-(28) 
by applying equations (25)-(28) to each family of stiffeners with the attributes given in Table 9. 
In the expressions that follow, the stiffener extensional modulus, shear modulus, eccentricity, 
moment of inertia, torsion constant, and transverse-shear correction factors refer to the 
corresponding effective quantities defined in Appendix A for a nonhomogeneous, specially 
orthotropic beam. The subscript and superscript "s" refers to the stringers and the subscripts and 
superscripts "dl" and "d2" refer to the two nonidentical families of diagonals. 
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D 26 = D 26 + H 3 E dl I dl - E d2 I d2 + i(G dl J dl - G d2 J d2 ) 2 

2L - v ' 4H 


plate G S J S BH 

D 66 = D 66 + — — + — j E dl I dl + E d2 I d2 + 

4H 4L 


(G dl J dl + G d2 J d2 j 
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4H B 
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A 45 — A 45 + — ^E dl A d |X z — E d2 A d2 x z J 


plate 4 E s A s x z L + (E dl A dl Xz + E d2 A d2 x z 2 )B 

A = A cc + 

55 55 4HL 


where the inplane-shear-deformation parameters are defined as 


s k Y G s 

Tv = .. 


and the transverse-shear-deformation parameters are defined as 


s k z G s 

X, = — ^ 

E s 


di _ k z G dl 


d2 _ k z G d2 


Setting these shear-deformation parameters equal to zero eliminates the effects of stiffener shear 
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deformation. 

/o' 

For the special case of equilateral- star cells with identical diagonals B = L and H = v 9 L 
and the stiffnesses reduce to 


plate 2 V3 E s A 
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^3E d A d / 
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A 22 = A 22 + 


4L 


I t m 
1 +: Ty 
3 


A - A plate 

AY 26 “ AV 26 


A - A plate ^sAJy ^AEjAj/ Ty^ 

A66 “ Ae6 + 6L + 4L [ Tj 


(G29) 

(G30) 

(G31) 

(G32) 

(G33) 


plate 2VTE s A s z s v^E d A d /_ d= \ 

B 1 1 — B , | + ,,, + — (z d + 3x Y z d | 


3L 


12L 


plate 73E d A d / d = 

B 12 — B 12 + (z d X Y Z d 


B i6 — B 16 


B 22 “ B 22 + 


plate . 3V3E d A d |_ . T v z 


4L 


Za + 


d = 
Y^d 


_ -r»P late 

^26 — E>26 


6L 
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_ plate 2V3'E S I S 

11-11 3L 12L 


(E d I d + 3G d J d ) 
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d 16 = d7 


D 22 =D r + ^E d i d+ G ^ 
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A 55 = Af + ^(2E s A s x a + E d A d x‘) 


where the subscript and superscript 11 d" indicates the properties of the diagonal stiffeners. 
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Appendix H 

Equivalent-Plate Stiffnesses for a Hexagon-Cell-Core Sandwich Plate 


The equivalent-plate stiffnesses for a sandwich plate composed of two different, anisotropic, 
laminated-composite face plates and the hexagon-cell core, made of a homogeneous orthotropic 
material, shown in figure 25 are presented in this appendix. The thicknesses of the bottom face 
plate (plate no. 1), the top face plate (plate no. 2), and the core are denoted by h p h„ and h c , 
respectively. The plate reference plane is located at the midplane of the core and the eccentricities 

of face plates no. 1 and no. 2 are given by e, = - i(h , + h c ) and e 2 = i(h 2 + h c ) , respectively. 

The core is composed of two member types. One member type is aligned with the y-axis of the 
plate, as shown in figures 21 and 25, and is referred to herein as a rib. The other member type 
makes an angle O with the x-axis and is referred to herein as a diagonal. Member attributes 
associated with the ribs and diagonals are indicated subsequently with the subscripts or 
superscripts "r" and "d," respectively. 

For a hexagon-core stiffener arrangement with each member made of a homogeneous 
orthotropic material, the effective moduli given by equations (A13a)-(A13c) reduce to the 
corresponding principal moduli of the orthotropic material, where the principal orthotropic- 
material axes are coincident with the X, Y, and Z beam coordinate axes shown in figure 6. Thus, 
the symbols E r and E d denote the extensional modulus E x and the symbols G r and G d denote 

the mean principal shear moduli ^G XZ G X ^' of the ribs and diagonals, respectively. Similarly, the 
stiffness-weighted eccentricities z r , z d , % , and z d vanish for homogeneous stiffeners centered on 

the plate reference plane. Moreover, the stiffness-weighted moments of inertia reduce to the 
corresponding second moments of area and the stiffness-weighted product of inertia vanishes for 
the rectangular stiffener cross sections. The cross-sectional areas of the core members are given 
by A r = h.t, and A d = h.t d , where t r is the rib thickness and t d is the diagonal thickness. The 

3 

moments of inertia associated with out-of-plane bending are given by I r = for the ribs and 

I d = for the diagonals. The planar geometric dimensions of the core, as shown in figure 21 

and 25, are given by the symbols a, b, c, and L. The symbols G J, and G d J d are the torsional 
stiffnesses of the ribs and diagonals, respectively. For a rib with a rectangular cross-section made 
of a homogeneous orthotropic material, equation (A 16) yields 


Gdr = G 
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fh c 
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(H2) 


J XZ J 


Likewise, the symbols k z and k d z denote the stiffness-weighted transverse-shear correction 
factor k z A /^' X7 of the ribs and diagonals, respectively, and k r Y and k d Y denote the stiffness- 

weighted inplane-shear correction factor k Y a /° xy of the ribs and diagonals, respectively. 

V G XZ 

For the notation just given, the equivalent-plate-stiffness contributions of the core to equations 
(47) are given by 


A n 


E d h c t d a 3 
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E d h c t d G d J d b 

19 + 2 

iz a 


(H9) 


j^core _ a b 

U 12 — , >. 3 

fb + c)L 


E h\ 

J -'d 11 c L d t 

'" J d J d 


(H10) 


136 



core _ 

D 16 =0 


(HI 1) 


. core C 

D 22 = 


r(b 


able 


E,h 3 c t r , b 4 I E d h 4 t d , ^ T a 2 ^ 

9 A ‘ 3 19 ' u il j ti 2 

24 CL \ 12 by 


core _ 

D 26 =0 


(HI 2) 
(HI 3) 


D 


core 

66 



E d hct d 

12 


+ 


G,J,L c 
8a 2 b 2 



(H14) 


. core H. 2 r 

A 44 = , . 2E d t d x z b + E r t r x z cL 

2a(b + c)LL J 

(H15) 

. core ^ 

a 45 =0 

(HI 6) 

core E d h c t d X z a 

55 " (b + c)L 

(HI 7) 

where the inplane-shear-deformation parameters are defined as 


r k Y G r 

T< - E, 

(HI 8) 

d k Y G d 

X Y - 

(HI 9) 

and the transverse-shear-deformation parameters are defined as 


r k z G r 

Xy = — ^ 

E r 

(H20) 

_ k z G d 
z ~ E 

J-'d 

(H21) 

Setting these shear-deformation parameters equal to zero eliminates the effects of stiffener shear 
deformation. The overall stiffnesses of the sandwich plate are obtained by substituting equations 


(H3)-(H17) into equations (47). 
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For the special case of a homogeneous, isotropic core comprised of hexagons with all sides 
having length L, the angle <F shown in figure 21b is 30 degrees, a = , b = L/2, and c = L. 

In addition, E =E , = E and G = G=G, where E and G are the extensional and shear moduli 
of an isotropic material, and k Y = k d Y = k' y = k z = k . For this case, equations (HI) and (H2) are 
expressed as 

Gd r = ^(l-K r ) (H22) 


G d J d = 


Gtjh c 


( 1 - K d) 
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The equivalent-plate core stiffnesses reduce to 
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Appendix I 

Equivalent-Plate Stiffnesses for an Orthogrid-Core Sandwich Plate 


The equivalent-plate stiffnesses for a sandwich plate composed of two different, anisotropic, 
laminated-composite face plates and the orthogrid core with homogeneous, orthotropic members 
shown in figure 26 are presented in this appendix. The thicknesses of the bottom face plate (plate 
no. 1), the top face plate (plate no. 2), and the core are denoted by h p h 2 , and h c , respectively. The 
plate reference plane is located at the midplane of the core and the eccentricities of face plates 

no. 1 and no. 2 are given by e, = - y(h , + h c ) and e 2 = i(h 2 + h c ) , respectively. The stringer 

and rib elements of the core are aligned with the x- and y-axes of the plate shown in figure 26, 
respectively. Member attributes associated with the ribs and stringers are indicated subsequently 
with the subscripts or superscripts "r" and "s," respectively. 


For an orthogrid core with each member made of a homogeneous orthotropic material, the 
effective moduli given by equations (A13a)-(A13c) reduce to the corresponding principal moduli 
of the orthotropic material, where the principal orthotropic-material axes are coincident with the 
X, Y, and Z beam coordinate axes shown in figure 6. Thus, the symbols E r and E s denote the 
extensional modulus E x and the symbols G r and G s denote the mean principal shear moduli 
\/G xz G xv °f the ribs and stringers, respectively. Similarly, the stiffness-weighted first moments 


of area z r , z s , z r , and z s vanish for homogeneous stiffeners centered on the plate reference plane. 


Moreover, the stiffness-weighted moments of inertia reduce to the corresponding second 
moments of area and the stiffness-weighted product of inertia vanishes for the rectangular 
stiffener cross sections. The cross-sectional areas of the stringers and ribs are given by A s = h c t s 
and A r = h c t r , respectively, where t s is the stringer thickness and t. is the rib thickness. The 


moments of inertia associated with out-of-plane bending are given by i s 



for the stringers 


and i r = | A f° r the ribs. The symbols L y and L x denote the spacing of the stringers and ribs, 

respectively, as shown in figure 26. The symbols G s J s and G r J r are the torsional stiffnesses of the 
stringers and ribs, respectively. For a rib with a rectangular cross-section made of a homogeneous 
orthotropic material, equation (A 16) yields 
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xz 3 


96 t s I G 


’h„ G 


2 


p= 1, 2, 3, 


P 


f 

tanh 

l 
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( 12 ) 


Likewise, the symbols k' 7 and k 7 denote the stiffness-weighted transverse-shear correction 
factor k z A A' X7 of the stringers and ribs, respectively, and k s Y and k y denote the stiffness- 

weighted inplane-shear correction factor k Y , A( XY of the stringers and ribs, respectively. 

V ^xz 

Based on the notation just given, the equivalent-plate-stiffness contributions of the core are 
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(13) 


(14) 


(15) 


As a special case, consider a sandwich plate with face sheets and a core are made from three 
different homogeneous, specially orthotropic materials. The core is made from a monolithic piece 
of orthotropic material. The face plates are relatively thin and, as a result, the transverse-shearing 
deformations of the face plates are presumed to be negligible. For this case, the principal material 
coordinate frames are taken to be identical to the structural coordinate frame shown in figure 26. 
The orthotropic-material properties are denoted by the extensional moduli E x and E y ; the shear 
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moduli G x , G xz , and G z ; and the Poisson’s ratios v x , and v x . In addition, superscripts are applied 
to these material-property designations to indicate their association with a given face plate or the 

core. For the stringer elements of the core, which are aligned with the x-axis; E s -* E x "' e , 
k Y G 5 -* k“ "G x ° re , and k z G s -* k)‘ r ’G“ re . For the rib elements of the core, which are aligned with 

the y-axis; E r E y , k Y G r -* k x G xy , and k z G r -> k 2 G yz . For a stiffener with a rectangular 
cross-section made of a homogeneous orthotropic material, equations (II) and (12) give 
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Using this information with equations (13) - (15) and (47) gives 
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Appendix J 

Equivalent-Plate Stiffnesses for a Star-Cell-Core Sandwich Plate 


The equivalent-plate stiffnesses for a sandwich plate composed of two nonidentical, 
anisotropic, laminated-composite face plates and the isosceles-star-cell core made of a 
homogeneous orthotropic material shown in figure 27 are presented in this appendix. The 
thicknesses of the bottom face plate (plate no. 1), the top face plate (plate no. 2), and the core are 
denoted by h p h 2 , and h c , respectively. The plate reference plane is located at the midplane of the 

core and the eccentricities of face plates no. 1 and no. 2 are given by e, = - J, (h , + h c ) and 

e 2 = i(h 2 + h c ) , respectively. The core is composed of two member types. One member type is 

aligned with the x-axis of the plate, as shown in figures 23 and 27, and is referred to herein as a 
stringer member. The other member type makes an angle with the x-axis and is referred to 
herein as a diagonal. Member attributes associated with the stringers and diagonals are indicated 
subsequently with the subscripts or superscripts "s" and "d," respectively. 

For an isosceles-star-cell core with each member made of a homogeneous orthotropic material, 
the effective moduli given by equations (A13a)-(A13c) reduce to the corresponding principal 
moduli of the orthotropic material, where the principal orthotropic-material axes are coincident 
with the X, Y, and Z beam coordinate axes shown in figure 6. Thus, the symbols E s and E d 
denote the extensional modulus E x and the symbols G s and G d denote the mean principal shear 
moduli N ( ; ( ; of the stringers and diagonals, respectively. Similarly, the stiffness-weighted 

first moments of area z s , z d , g s , and i d vanish for homogeneous stiffeners centered on the plate 

reference plane. Moreover, the stiffness-weighted moments of inertia reduce to the corresponding 
second moments of area and the stiffness-weighted product of inertia vanishes for the rectangular 
stiffener cross sections. The cross-sectional areas of the core members are given by A s = h c t s and 
A d = h.t d , where t s is the stringer thickness and t d is the diagonal thickness. The moments of 

3 3 

inertia associated with out-of-plane bending are given by I s = jy f° r the stringers and l (l = 

for the diagonals. The planar geometric dimensions of the core, as shown in figure 23 and 27, are 
given by the symbols B and H, where B and H are the base and height, respectively, of either 
isosceles triangle forming the star shape. The symbols G J s and G d J d are the torsional stiffnesses 
of the stringers and diagonals, respectively. For a stringer with a rectangular cross-section made 
of a homogeneous orthotropic material, equation (A 16) yields 
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Likewise, the symbols k z and k z denote the stiffness-weighted transverse-shear correction 
factor k z * /fps of the stringers and diagonals, respectively, and k Y and k ' denote the 

stiffness-weighted inplane-shear correction factor k Y A A' xv of the stringers and diagonals, 
respectively. 

Based on the notation just given, the equivalent-plate-stiffnesses contributions of the core are 
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where the inplane-shear-deformation parameters are defined as 

s k Y G s 

T -—r 

_ k Y G d 
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J-'d 

and the transverse-shear-deformation parameters are defined as 



(J 18) 
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The overall stiffnesses of the sandwich plate are obtained by substituting equations (J3)-(J 1 7) into 
equations (47). 

For the special case of an isotropic star-cell core comprised of triangles with all sides having 
length L and thicknesses t s = t d = t, the angle O shown in figure 23b is 60 degrees, H = ,/G L , and 
B = L. In addition E s = E d = E and G = G d = G, where E and G are the extensional and shear 
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moduli of an isotropic material, and k Y = k d Y = k s z = k z = k . For this case, equations (Jl) and (J2) 
are expressed as 


G S J S = G d J d = 


GtV 


(l-K) 


( 122 ) 
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Thus, the equivalent-plate core-stiffness contributions become 
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